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ABSTRACT 

This  study  treats  the  effect  of  the  interaction  between  underground 
structures  and  the  surrounding  soil  in  reducing  the  loads  transmitted  to  the 
structure,  the  so-called  "arching"  phenomenon. 

A  continuum  theory  of  soils  proposed  by  G.  A.  Geniev  is  applied  to  a 
quasi-static,  plane-strain  problem  of  arching.  The  basic  partial  differential 
equations  are  shovn  to  form  a  hyperbolic  set  end  are  solved  by  the  method  of 
characteristics.  Consistent  stress  and  velocity  fields  are  obtained. 

Comparison  with  available  experimental  results  shows  that  the  Geniev 
theory  underestimates  the  surface  pressure  required  for  failure  of  an 
underground  structure  in  relatively  dense  granular  soils.  The  source  of 
this  difficulty  is  explained  and  an  aonrovlmate  method  of  overcoming  it  is 
presented. 

A  simplified  extension  to  a  theory  taking  account  of  Inertia  of  the  soil 
and  unsteady  motions  is  treated  in  an  appendix. 
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I .  INTRODUCTION  TO  THE  SOIL- STRUCTURE  INTERACTION  PROILEH 
I . 1  Object  and  Scope 

The  purpose  of  the  main  part  of  this  report  (Chapters  1*6)  is  to 

investigate  the  applicability  of  a  continuum  theory  of  granular  media  to 

the  problem  of  soil-structure  interaction  of  underground  structures-  The 

I  2 

problem  essentially  involves  a  determination  of  the  effect  of  ••arching" 
in  granular  soils  on  the  pressure  transmitted  to  a  buried  structure  when  the 
ground  surface  is  subjected  to  a  loading  in  the  form  of  an  air  overpressure 
such  as  would  result  from  the  detonation  of  a  high-yield  nuclear  device- 

The  air  overpressure  loading  is  assumed  to  result  from  an  air  burst 
occurring  either  high  above  the  ground  surface  or  near  the  ground  surface 
with  a  euperseismie  velocity  of  propagation  of  the  pressure  wave,  so  that 
only  the  air-induced  pressure  nead  be  considered.  The  al r-over-pressure 
loading  is  assumed  to  vary  tlcatly  with  time  and  to  be  uniformly  distributed 
over  an  infinitely  large  area  as  far  as  the  underground  structure  is  concerned. 
A  pressure  distribution  of  this  type  is  adequate  for  an  introductory  study. 
for  the  truly  dynamic  problem  a  more  realistic  pressure  distribution  is 
needed.  In  the  superset amic  case  the  pressure  wave  propagates  doirwird 
through  the  soil  Inclined  at  an  angle  to  the  ground  surface;  3  harever  for 
simplicity  the  pressure  wave  is  assumed  to  propagate  vertically  bom  ward 
with  the  weve  front  parallel  to  the  gr«ntnd  surf-stc- 

The  soil  is  considered  to  be  a  homogeneous.  Isotropic,  granular 
material  having  some  cohesion  and  satisfying  the  Cmulomb-hebr  failure 
criterion  which  states  that  slip  or  failure  occurs  at  a  point  when  the  shear 
stress  on  we  place  through  the  point  is  equal  to  the  um  ef  the  cohesion 
and  a  function  of  the  normal  stress  acting  on  the  same  plane,  bn  no  plane 


I 


can  the  shear  stress  exceed  this  sum.  The  function  of  the  normal  stress  is 


taken  to  be  a  constant,  namely,  the  tangent  of  the  angle  of  internal  friction 
times  the  normal  stress.  The  fact  that  the  failure  criterion  depends  on 
the  hydrostatic  part  of  the  stress  tensor  and  not  simply  on  the  deviatoric 
part,  as  is  the  case  in  the  theory  of  plasticity  for  metals,  adds  an  Important 
complication  to  the  problem  as  Mill  be  seen  later. 

The  structure  to  be  considered  s  a  long,  rectangular  plate,  simply 
supported  at  the  edges.  It  is  assumed  that  the  length-width  ratio  of  the 
plate  is  large,  so  that  conditions  of  plane  strain  Mill  prevail  away  from 
the  ends  of  the  long  side.  It  follows  from  these  esstmptions  that  the  behavior 
of  the  structure  may  be  considered  to  be  that  of  a  beam,  and  the  problem  My 
be  formulated  in  two-dimensional  plane  strain. 

The  depth  of  burial  of  she  underground  structure  is  in  the  shallow 
range  with  the  meaimum  depth  of  Interest  on  the  order  of  twice  the  span 
length,  it  is  in  this  range  of  depths  that  the  effect  of  arching  In  the  soil 
is  expected  to  be  a  maxlmsm.  Although  the  attenuation  of  free-fleld  alr- 
ove'pressure  by  the  soil  My  be  significant  In  this  tenge.4  only  the 
redistribution  of  pressuiu  caused  by  s«i ’-structure  interact  ion  wi  I  <  be 
considered  in  determining  the  pressure  transmitted  to  the  structure- 

An  important  effect  of  the  soil  on  the  net ure  of  the  pressure  wove, 
particularly  for  dynamic  studies,  leads  to  en  increase  in  the  rise  iIm  of 
the  pressure  pulse  as  the  wave  propagates  downward  through  the  soil  5‘*.  This 
•nc reeved  rise  time  tends  to  reduce  dynamic  effects. 

1.2  Summary,  of  Previous  Work 

(Vie  Of  the  earliest  studies  of  the  soil-structure  interaction 
problem  to  appear  in  the  literature  wet  Tertaghi's  "Trap  boor"  experiment-17 
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With  h  simple  device  he  showed  that  a  small  downward  movement  of  an  under- 
ground  structure  results  in  a  large  decrease  in  the  intensity  of  pressure 
applied  to  the  structure.  This  type  of  behavior  is  of  interest  for  the  design 
of  tunnel  liners  or  large  underground  conduits  as  well  as  the  design  of 
underground  blast-resistant  structures. 

Several  methods  have  been  used  to  evaluate  the  effect  of  soil- 
structure  interaction  on  underground  structures  subjected  to  air-blast  load¬ 
ings.  One  study  was  made  by  considering  the  soil  to  be  an  ideal  elastic 
8 

material.  Spherical  and  cylindrical  models  were  treated  by  considering  a 
uniformly  distributed,  radially  symmetric  pressure,  and  using  equations  of 
equilibrium  and  compatibility  fron  the  theory  cf  elasticity.  The  method  was 
applied  to  shapes  of  structures  other  than  spherical  by  considering  the 
compressibility  of  the  structure,  that  is,  the  change  in  volume  due  to 
externally  applied  loads. 

Another  approach  considered  only  the  Rankine  passive  earth  pressure 

9 

developed  by  the  soil  as  a  flexible  two-hinged  arch  moved  Into  the  soil  mass. 

The  problem  was  reduced  to  a  single-degree-of-f reedom  system  consisting  o'  a 
lumped  mass  restrained  by  a  weightless  spring  with  the  configuration  being 
defined  by  the  radial  deflection  of  the  arch  at  the  quarter-point  rf  the  arch 
rib.  For  all  arches  considered,  the  response  was  found  to  occur  ii  the  plastic, 
or  yielded,  part  of  the  assumed  elast';-p*nstic  resistance  diagram  for  t'te 
so!  I . 

A  more  realistic  approach  utilising  the  classical  Rankine  earth 
pressure  theory  has  been  made  by  considering  both  the  active  and  passive  states 
of  stress  in  the  soil  and  the  different  soil  masses  entering  Into  the  response 
as  the  arch  deflects.3  In  this  cate,  the  structure  consisted  of  four  rigid, 
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massless  bars  connected  by  rotation  resisting  springs,  in  determining  the 

dynamic  response  of  the  two-degree-of-f reedom  system,  the  masses  were 

assumed  to  be  concentrated  at  the  bar  connections.  The  blast  overpressure 

acting  on  the  ground  surface  of  the  soil  failure  wedge  was  taken  into 

account  in  determining  the  total  resistance  of  the  soil  mass. 

A  dynamic  analysis  based  on  the  formation  of  vertical  slip  planes 

in  the  so>l  between  the  structure  and  the  ground  surface  and  a  uniform 

distribution  of  vertical  stress  on  horizontal  sections  has  been  made  for 

10 

pseudo-steady  state  air  overpressure.  In  this  analysis  the  failure  mass 
of  soil  is  considered  to  undergo  a  rigid  body  motion  along  with  the 
structure-  tn  this  theory  no  account  is  taken  of  the  propagation  of  the 
pressure  wave  through  the  soil  or  the  associated  dynamic  soil-structure 
interaction. 
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A  list  of  primary  symbols  used  in  the  first  six  chapters  is 

summarized  for  convenience  in  the  following. 

a  .a  ■  components  of  acceleration  in  the  x  and  y  directions 
a  y 

H  *  depth  of  burial  of  structure 

K  «  ratio  of  horizontal  to  vertical  stress  at  a  point  in  the  soil 
k  ■  soil  cohesion  in  units  of  stress 

L  -  span  length  of  structure 

pQ  *  intensity  of  overpressure  at  the  ground  surface 

q  *  intensity  of  pressure  loading  on  the  structure 

r,6  *  polar  coordinates 

v  ,v  *  components  of  velocity  in  the  x  and  y  directions 
x  y 

X,Y  "  intensities  of  body  force  in  the  x  and  y  direction 
x,y  ■  rectangular  coordinates 

3  ■  direction  of  the  velocity  vector  measured  counterclockwise 
from  the  positive  x-exis  in  rectangular  coordinates 

a  m  direction  of  the  velocity  vector  measured  counterclockwise 
from  the  positive  radius  vector  in  polar  coordinates 

0  ■  direction  of  the  maximum  principal  stress  measured  counter* 

clockwise  from  the  positive  x-exis  in  rectangular  coordinates 

$  ■  direction  of  the  maximum  principal  stress  measured  counter¬ 
clockwise  from  the  positive  radius  vector  in  polar  coordinates 

r  »  */4  -  cp/2.  Physically,  the  Inclination  of  possible  slip 
planes  to  the  direction  or  princloa’  »t>uss 

c  ,c  -  components  of  strain  in  the  x  and  y  directions 
x  y 

i  ,t  •  components  of  strain-  rete  in  the  x  and  y  directions 
x  y 

T  ■  shear  strain 

xy 

r  »  shear  strain  rate 

ny 
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1,1}  ■  rectangular  coordinates  on  the  characteristic  plane 

p  «  density  of  the  soil  mass  per  unit  volume 

o  ,o  •  components  of  stress  in  the  x  and  y  directions 

x  y 

x  ■  shear  stress 

xy 

Oj,o2  «  principal  stress 

»  normal  stress  and  shear  stress  on  failure  plane 

rs  n 

civ  «  combined  stress  and  velocity  variables 

cp  »  angle  of  internal  friction  of  the  soil 
X  *  angle  between  the  x-axis  and  the  charccterl st! c  curve 


2.  DYNAMIC  THEORY  OF  GRANULAR  MEDIA 


2. !  Assumptions 

The  theory  of  the  dynamics  of  granular  media  to  be  presented  in 
this  and  the  next  two  chapters  follows  that  of  G.  A.  Geniev".  In  applying 
this  theory  to  the  problem  of  soil  dynamics,  it  is  assumed  that  the  soil  will 
benave  as  a  continuous  and  homogeneous  granular  medium  having  a  shearing 
resistance  consisting  of  the  cohesion,  which  is  assumed  to  be  a  constant 
independent  of  the  stress  state  at  a  point,  and  internal  friction,  which  is 
assumed  to  be  a  linear  function  of  the  mean  normal  stress. 

It  has  been  shown  by  experiment1  that  the  resistance  to  slip  along 
a  plane  through  a  typical  point  in  the  medium  is  given  by 

It  I  ■  o  tan  a  +  k  (2.1.1) 

n  n 

where  o  , t  “  the  normal  and  tangential  components  of  the  stress  vector 
n  n 

acting  at  the  point, 

®  ■  the  angle  of  internal  friction, 
k  ■  the  cohesion. 

Planes  in  the  soil  medium  on  which  the  relation 

max  £|t|  -  (ofl  tan  *  +  k)  j  ■  0  (2.12) 

1s  satisfied  are  called  slip  planes.  At  any  point  In  a  region  of  slip  there 
are  two  planes  along  which  slip  may  occur.  The  one  along  which  slips  actually 
occur  is  called  the  "active  slip  plane." 

When  egul librium  is  destroyed,  a  velocity  field  resulting  from  the 
relative  motion  of  the  toil  particles  it  set  up  In  the  toil.  It  Is  assumed 
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that  this  velocity  field  will  satisfy  the  equation  of  continuity  from  the 

theory  of  continuous  media.  It  is  further  assumed  that  when  How  occurs 

the  stress  components  may  be  expressed  in  terms  of  the  strain  rates. 

As  stated  in  Chapter  1,  the  two-dimensional  plane  problem  will 

be  considered  in  which  the  stress  components  a  ,  a  ,  x  and  the  velocity 

x  y  xy 

components  v  and  v  are  expressed  as  functions  of  the  space  variables  x  and 
x  y 

y  and  the  time  t-  Since  relatively  large  displacements  are  expected,  the 
problem  has  been  formulated  in  Eulerian  coordinates  so  that  the  motion  is 
expressed  in  terms  of  the  components  of  the  velocity  field  rather  than  the 
displacements  of  the  particles.  After  a  complete  velocity  solution  has  been 
obtained,  the  particle  displacements  may  be  found  by  direct  integration. 

2.2  Equations  of  Hot  ion 

Figure  2.1  shows  an  infinitesimal  elemant  of  the  soil  modi  urn  which 

is  in  a  state  c#  flow  and  which  it  acted  upon  by  the  direct  stresses  a x  and 

?  and  the  shear  stress  t  The  sign  convention  adopted  here  takes  com- 
V  *v 

pressive  stresses  positive,  so  that  the  inward  normal  stress  acting  on  a 
surface  of  the  element  shown  in  Fig.  2.1  is  positive  if  the  inward  normal  to 
the  surface  is  in  the  direction  of  the  negative  coordinate  axis.  Positive 
shearing  stress  on  a  surface  acts  In  the  negative  direction  of  one  coore  nate 
axis  when  the  Inward  normal  to  the  surface  acts  in  the  negative  dlrectlto 
of  the  other  coordinate  axis-  This  sign  co.-.vt.r,!on  <*  opposite  to  the  usual 
convention  In  the  theory  of  elasticity;  hewever,  it  is  convenient  since 
the  normal  stresses  usually  encountered  in  soils  are  compressions. 

If  X  and  V  are  the  intensities  of  the  body  force  and  p  Is  the  mass 
per  unit  vo  I  we ,  the  equation  of  motion  in  the  x-dl  rest  Ion  of  the  element  In 
Fig.  1.1  is: 


0  X  dx  dy  -  (ox  +  37*  dx)  dy  +  ox  dy  -  (txy  +  dy)  dx 

+  t  dx  •  p  dx  dy  a 
xy  w  '  x 

where  ax  is  the  acceleration  of  the  element  in  the  x-di rection.  This  equation 
may  be  written  at 

Similarly,  the  equation  of  motion  in  the  y-di rection  yields 


In  Culerian  coordinates  the  acceleration  terms  a  and  a  in  Cqs-  (2-2.1)  and 

x  V 

(2.2.2)  are  expressed  as  functions  of  the  velocity  components,  v  and  v  by 

h  y 


dvB 

dv_ 

dv_ 

ST* 

V*ST  + vy3T 

(2.2.3) 

dv 

dv 

3t*  + 

y«3T  +  W 

(2.2.4) 

Using  these  expressions,  the  equations  of  equilibrium  bi 


do  dr  dv  dv  dv 

I  ,  x  .  _  x  .  ..  x  .  .. 


* ;  (3T  *  V  "St®*  v*5£  +  \3T  (,  ,  5) 


,  dr  do  dv  dv  dv 


C.2.«) 


J  J  N»i.urt  Criterion 

At  stated  earlier,  slip  will  occur  when  the  stress  components 
satisfy  the  relation 


m»x  |  |tJ  -  (o^  tan  9  +  k)  |  -  0 


The  points  on  a  Mohr  stress  diagram  which  satisfy  this  relation  form  an 
envelope  of  possible  states  of  stress  at  failure.  One  such  stress  state 
is  shown  on  the  Coulomb-Mohr  diagram  In  Fig.  2.2.  The  radius  of  a  Mohr 
circle  for  a  state  of  stress  satisfying  the  failure  criterion  is 
4  (a,  -  a,)  ■  t  where  o.  and  a.  are  the  maximum  and  minimum  principal 
stresses.  The  principal  stresses  a(  and  Oj.  and  t  are  related  to 
xy- components  of  stress  by 

O,  -  {  (ox  ♦  «y>  ♦  {  <2-3- » 

°2-i(ox  +  V  'I  J(ax‘  V*  +4T*yl 

Tma».  -  I  iol  *  *2>  “  \  *  4 

expressions  for  e  and  t  are  obtained  from  Fig.  2.2  to  be 

ft  ft 

°n  "  1  ion  *  V  *  I  >/(,«  •  ^  *'"•  •  (2*'*4) 

\mi  -  •/ 44  V  “*•  •  (,  ,*>) 

Substituting  Cos.  (23.4)  and  (2.3.1)  Into  tbe  Coulomb  failure 
criterion,  tg.  (2.1.1),  one  obtains: 

(*M  *  Oy)2  ■<  4  t^2  •  s»n2f  (u^  ♦  Oy  ♦  2k  cot  f)2  (2.3.C 
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(2.3.2) 

(2.3.3) 


It  may  be  seen  from  Fig.  2.2  that  the  two  possible  slip  planes  are  inclined 
at  angles  ±7  to  the  direction  of  the  maxim'*  principal  stress,  where  the 
angle  7  is  given  by 


7  f  (2-3.7) 

The  angle  between  the  slip  planes  and  the  x-axis  is  determined  by 
introducing  the  angle  0  which  is  the  angle  between  the  direction  of  maximum 
principal  stress  and  the  x-axis.  Then  the  slip  planes  will  form  at  angles 
0  t  7  to  the  x-axis,  as  is  shown  in  Fig.  2.3. 

2.4  Equation  of  Continuity 

According  to  the  theory  of  continuous  media,  the  velocity  components 
and  the  density  are  related  by  the  principle  of  conservation  of  mass  In  the 
fol  lowing  way: 

|f  the  medium  is  incompressible,  the  equation  of  continuity 

becomes: 


(2.3») 


2.S  The  Kinematic  Mining 

The  kinwMtlc  relations  are  expressions  which  relate  the  stresses 
to  the  velocity  field.  The  two  such  relations  to  be  considered  here  a»e 
based  on  the  following  assumptions: 

I.  The  sheer  strain  rate  Is  taken  to  be  a  maxims*  along  the 
direction  of  active  slip,  (2.1.1) 

II 


2.  Tne  velocity  vector  is  restricted  to  coincide  with  the 


directions  of  active  slip,  (2.S.2). 

The  first  of  these  relations  considers  the  soil  to  be  essentially 
a  truly  continuous  medium  while  the  second  attempts  to  account  for  the 
behavior  of  a  medium  consisting  of  a  system  of  particles  whose  dimensions 
are  small  but  not  insignificant- 

The  theories  resulting  from  each  of  tliete  relations  will  he 
referred  to  as  the  continuum  theory  and,  as  Geniev  calls  it,  the  macro- 
structural  theory,  respectively. 


2-5-1  Kinematic  Relation  for  the  Continuum  Theory.  Figure  2.4 
shows  the  Mohr  circle  of  the  strain  rates  4  ,  <  ,  and  -j  7  which  are 

*  y  *  *7 

related  to  the  velocity  coaponents  by 


and 


The  direction  of  the  maaimym  shear  strain  rate  with  respect  to  the 
x-axis  is  given  by  the  angle  *,  which  from  Flo-  2.4  is 


tan  2* 


(2*1' 
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Since  the  maximum  shear  strain  rate  occurs  along  a  slip  line  in 


this  theory,  the  angle  X  is; 


X  -  (p  i  y) 

In  Fig.  2.5,  this  relation«hip  is  shown  clearly.  Then  from  Fig.  2. 


tan  2(P  t  y) 


From  the  Mohr  circle  of  stresses  shown  in  Fig.  2.5, 


Combining  tqs-  (2. 5. 2)  and  (2.5. 3) ,  the  kinematic  relation  between 
and  velocity  fields  it  obtained  for  the  continuum  theory. 


dv_ 


This  may  also  'oe  written  at  foli»*s: 
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(2-5.2) 

(2.5.3) 

the  stress 

(25.4) 

(2.1.5) 


2.5.2  Kinematic  Relation  for  the  Macro  Structural  Theory.  The 
second  kinematic  relation  is  based  on  the  assumption  thet  the  velocity  vector 
coincides  with  the  direction  of  active  slip.  This  assumption  follows  by 
considering  a  state  of  motion  in  which  the  soil  grains  are  moving  In  tubes 
or  layers  of  flow  formed  by  adjacent  stream  lines  as  shown  in  Fig.  2.7.  The 
smallest  spacing  of  two  stream  lines  is  given  by  the  dimensions  of  e  single 
particle,  and  the  velocity  of  flow  in  any  layer  varies  with  the  width  of  that 
layer.  If  relative  motion  occurs  in  the  soli  medium,  slip  will  take  place 
along  the  boundaries  of  the  layers  of  flow,  that  is,  along  the  stream  lines. 
Since  the  tangent  to  a  stream  line  represents  the  direction  of  the  velocity 
vector,  the  velocity  vector  may  be  assumed  to  coincide  with  the  direction 
of  lines  o*  active  slip. 

In  Fig.  2.7  a  field  of  active  slip  lines  and  the  direction  of  the 
velocity  vector  at  a  typical  point  is  shorn, 

If  o  Is  the  angle  between  the  direction  of  the  velocity  and  the 
a-eais,  Fig.  2.7  gives 

v 

tan  O  •  -JL  (2.S.4) 

a 

Here,  the  direction  of  the  slip  line  at  the  point  Is  also  given  by  the  angle  o, 

«  •  0  *  t  {J.4.J) 

Then 

ten(0*T).^£  (l$.«> 

In  lg.  (2.5.J),  the  stress  components  ere  related  to  the  angle  20  the 
Hohr  circle. 
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28  -  -- 

°Y 


tan  2(8  t  r>  -  -  2  t<n  iP  *  Xl  _ _ !*__ 

I  -  tan2  (8  t  r)  ,  . 


(2.5.9) 


tan  2(8  *  r)  -  -r^"  *  V**  2 
I  *  tan  28  tan 


(2.5.10) 


an  expression  for  tan  28  My  be  obtained. 


ton  18  * 


,  *  »  , 

I'r-r1  T3T57 


V  ¥ 


l  1 1  ( Jt  .  JS)  - ' 

*  F 


(1.5.11) 


Cooiinin«  t*«.  (I.s.ll)  and  (2.5.2),  and  notlitf  that  ana  obtain*  an  <4»*tlan 
wbisti  ratata*  tba  »tr**»  and  velocity  field* 


,  i  £  *  'r*  * s"  • 

ji _ *  » 


I.  -  *  .  *  y 

'  »I  <»* 

®  T 


(1.5.12) 


Thus,  the  five  unknown  functions  ox>  *Xy«  vx  anc*  vy  ar*  9'ven 
by  the  f :  ve  equations.  (2.2.5),  (2-2.6).  (2.3.6).  (2.3-8)  or  (2.3.9).  and 
(2.5.5)  or  (2.5  12)  . 

!t  may  be  shown  that  the  kinematic  relations  coincide  in  any 
region  where  the  streamlines  are  parallel  and  the  flew  is  steady.  In  this 
case  the  velocity  vector  does  not  vary  with  time  or  position,  that  is. 

-w-5dt+Sd**^dy-0. 


where  dx  and  dy  refe--  to  a  segment  of  streamline. 

Since  the  flow  is  steady.  -  0.  Oividing  the  above  equation  by  dt.  on. 
obtains  the  fol loving  expression. 


Vy  -  0 


(2.513) 


The  velocity  v  may  be  expressed  in  terms  of  v^  and  vy  by: 


2  2  2 
V  •  V  ♦  * 

x  y 


(2  5.  U) 


Differentiating  lq.  (2  5.1*)  -*U  r*Wect  to  x. 
Obtained.  Similarly^-  is  found. 


tor  •» 


V-  t 

£  "  v  (y«  W  "y 
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(25.15) 


Substituting  Eqs.  (2.5.  IS)  into  (2.5. >9)  and  solving  the  resulting 
v  v 

quadratic  expressions  for  and  — ,  it  turns  out  that: 

v  v 

*  y 


dv  civ 

”» -  &L :  ?i 

•'  IT 


(2.5.16) 


Tnus,  it  may  be  seen  that  the  two  kinematic  relations  given  by  Eqs.  (2.5.5) 
and  (2.5.12)  coincide  in  any  region  of  steady  flow  with  parallel  stream  lines. 
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3.  TRANSFORMATION  OF  EQUATIONS  Bv  THE  METHOD  OF  CHARACTERISTICS 


3 . 1  Method  of  Solution 

Since  the  five  equations  relating  the  five  unknowns  form  a  system 
of  first  order,  quasi-linear  partial  differential  equations  of  the  hyperbolic 
type,  it  is  possible  to  formulate  a  procedure  for  the  solution  of  these 
equations  by  the  method  of  characteristics**'**’*^.  In  this  formulation,  the 
system  of  five  equations  in  five  unknown--  is  replaced  by  the  differential 
equations  of  the  characteristics  and  the  differential  equations  called 
‘‘characteristic  equations,"  to  be  satisfied  along  the  characte  istics.  The 
solution  obtained  by  integrating  the  cha. acter i st i c  equations  along  the 
characteristics  is  equ  valent  to  the  solution  of  the  original  system,  but  is 
generally  more  convenient  to  obtain  since  the  characteristic  differential 
equations  and  the  differential  equations  of  the  characteristics  are  ordinary 
rather  than  partial  as  in  the  original  problem. 

The  method  will  be  developed  for  the  quasi 'Static  case  of  an 

dv  dv 

incompressible  soil  medium  In  this  case  the  terms  x/dt  end  y/dt  art 
neglected,  and  p  is  taken  as  constant,  which  results  in  the  five  equations: 


» 


da 

i  , _ *  . 

;  (3T  * 


♦  V 


3T 


(3.1  I) 


y 


2 

p 


df  do 


dv ,  dvy 

v  ■»  v  — 

*  dx  >  jy 


(3.1.3) 


(o  -  a  )  *  a  4  t  *  ■  sin*  9  (0  a  0  a  2k  cot  9)  *  (3.1.3) 

x  y  xy  w  x  y 


'ST 


Sr1 


(3.1.4) 


II 


tan  <p 


(3.1.5) 


2  LUL 


°x  '  0 


r  2  ay  o* 

T  m'& - — SC - Sv 


■)  tan  9 


L*i  . 
a  -  a 
*  y 


1  ^  *  T*  1  t-*  9 
-  x  I _ 

'  + 1  ^  Un  9 


(3.1.6) 


The  last  equation,  (3.1.5)  or  (3.1.6),  is  taken  depending  on  which  kinematic 
relation  is  used. 


3.2  Transformation  of  Variables 

It  is  convenient  to  define  four  new  unknowns  in  term*  of  the 

original  five,  to  that  the  algebraic  equation  of  the  failure  criterion  (3.1.3), 

is  automatical ly  satisfied.  This  will  reduce  the  original  system  of  five 

equations  to  four  equations  in  the  four  new  unknowns. 

The  principal  stresses  0|  and  0j  My  be  expressed  In  terms  of  the 

stress  components  o  .  u  ,  and  T  from  the  Nohr  circle  of  stresses,  Fig.  2.5: 
a  y  ay 


"i  •  s  (°.  *  v  *  t  ■/<».•«/**  >.7 

«!  •  T  ♦  V  •  W'V 


so  that 


and 


o  +o 

X  y 


o,  +  r2 


(3.2.2) 


Taking  the  square  root  of  both  sides  of  Eq.  (3.1.3)  and  using  Eqs.  (3.2.1) 
and  (3.2.2),  one  obtains  the  equation  of  the  failure  criterion  in  terms  of 
the  principal  stresses: 

a,  +  o,  -  -  2k  cot  (3.2.3) 
1  2  Sin  <p 


The  stress  components  o  ,  c  ,  and  t  are 

Ay  Ay 


«x  •  j  (°|  +  °2)  +  1 


\  *  2  (°l  *  °2)  *  1 


T«y  "  2  (°1  *  ,1) 


given  by: 


(Cj  -  o?)  cos  20 

(o,  -  o2)  cos  20 

20 


Using  Eq  (3.2.3),  these  equation*  become; 

(3.2  4) 

(3.2.5) 

(3.2.6) 

where 


-  oo  (»♦  sin  9  cos  20)  -  k  cot  9 

■  oo  (i  •  sin  9  cos  20)  -  k  co*  v 

t  •  po  (sin  9  sin  20) 
xy 


m 


On  the  Coulomb-Mohr  diagram,  po  represents  the  distance  from  the  intercept  of 

the  failure  envelope  and  the  normal  stress  axis  to  the  center  of  the  Hohr 

circle  as  shown  in  Fig.  2.2. 

The  velocity  components  v  and  v  may  be  expressed  as: 

x  y 

v^  ■  v  cos  a  (3 .2. 7) 

vy  ■  v  sin  a  (3.2.9) 

where 

a  ■  angle  between  the  x-exis  and  the  velocity  vector  v  (Fig.  2.6) . 

3 . 3  Characteristics  of  the  Continuum  Theory 

Using  Eqs.  (3.2-4).  (3.2.5),  (3.2  6).  (3.2.7)  and  (3.2.6),  the 
original  equations  with  the  kinematic  relation  of  the  continuum  theory  My 
be  expressed  in  terms  of  the  new  unknowns  e,  0.  o>,  and  v. 


^2  cos  (9  ♦  r)  ♦  sl«  (P  ♦  r>  ♦  2o  ten  f  ^  cos  (P  ♦  r)  ♦  sin  (P  ♦  r)  J  ♦ 

*  TSTi  [J  tc*  *  T)  *  S  (p  *  T) ]  *  • 0  <*•*•  '> 

cos  (P  -  r)  ♦  ^  sin  (p  -  r)  -  Jo  •"  r  cos  (p  -  r)  ♦  sin  (P  -  r)  j  * 

*  tsr;  [J  «•* (*  ♦  T)  *  S  * r)] 4  " 0  (3  ,  « 
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v  tan  j  Q  -  o  ±  r)  j  Fg  cos  (p  i  r)  +  g  sin  (P  ±  r) 


[g  cos  (P  t  r)  +  g  sin  (p  *  r) 


(3.3.3) 


v  cot 


a  -  (P  t  r) 


g  sin  (p  1  r)  -  g  eos 


(P  1  r) 


g  sin  (p  1  r)  -  g  cos  (P  t  r)J  -  o 


(3.3.4) 


where  A  ■  *  sin  (P  -  y)  *  Y  cos  (P  -  r) 

#  sin  (P  +  y)  +  Y  COS  (P  +  y) 

The  first  two  equations  of  this  set  represent  the  equations  of 
equilibrium  in  the  directions  perpendicular  to  the  slip  lines.  Cqs  (3.3.3) 
and  (3.3.4)  express  the  conditions  that  the  normal  strain  rates  vanish  In 
the  directions  along  and  perpendicular  to  the  slip  lines-  It  Is  now  easy  to 
see  that  the  set  (3.3.1),  (3.3.2)  is  equivalent  to  the  set  (3.1.1),  (3.1.2). 
Moreover,  (3.3.4)  and  (3.3.3)  taken  together  guarantee  that  the  voliaeetric 
strain  rate  it  tero  (*.he  center  of  the  circle  in  Fig.  (24)  Is  at  the  origi  i) . 
In  addition,  the  norma  I  strain  rate  of  tero  on  a  slip  line  then  lapllet  a 
maxi  mm  shear  strain  rate  on  this  line.  The  equivalence  *»l  th  the  orlglnai 
set  for  the  first  kinematic  condition  is  now  apparent- 

The  characteristics  of  Cqs-  (3.3.1)  to  (3.3.4)  nay  be  defined  as 
those  curves  y  »  y(x)  across  which  It  it  possible  for  the  derivatives  of  the 
unknown  functions  o,  P,  Q,  and  v  to  exhibit  finite  discontinuities  while  the 
functions  themselves  remain  continuous*3' **• 
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Along  the  curve  y  ■  y(x) ,  the  unknowns  satisfy  the  following 
differential  relations: 


do  '  15  dX  +  |y  dy  ;  dP  '  dx  +  ly  dy  ; 
^  "  lx  dx  +  ly  dy  ’  dv  "  |x  dx  +  ly  dy 


(3.3.5) 


Equations  (3.3.1)  to  (3.3.4)  and  (3.3.5)  are  a  system  of  tight 
algebraic  equations  in  the  eight  unknown  derivatives  of  o,  P,  a,  and  v  with 
respect  to  x  and  y.  Eliminating  all  derivatives  with  respect  to  x  from 
these  equations,  the  following  equations  are  obtained  Involving  only  derivatives 
with  respect  to  y: 


- -  — 

- 

•l 

b. 

cl 

0 

£ 

*. 

*2 

b2 

C2 

0 

• 

15 

• 

s 

0 

0 

c3 

d3 

£ 

s 

0 

0 

C4 

d4 

_v 

where 

a(  •  dx  sin  (0  ♦  y)  -  dy  cos  (P  ♦  ri 

Sj  «  dx  sin  (0  -  y)  -  dy  cos  (0  -  r) 

b|  »  2«  tin  f  •| 

b,  «  -2a  tan  g  a, 

*  1  * 


2) 


c3  -  V  tan 


COS  $  1 

a  ■  (M  r) 


c,  -  V  cot 
4 


dx 


d.  - 


j  Tdx  sin  (p  t  r)  -  «*y  «>*  (P  1  r)J 

-  (p  t  r)  j  |^dx  cos  (P  t  r)  +  dy  »'«  O  *  r)J 
sin  (Mr)  -  <*y  cos  (P  t  r)j 
n  (P  ±  r>] 


(3.3.7) 


dx  cos  (6  ;  r)  dy  sir 


cos  (p  +  r)  do  4  2o  tan  rp  cos  (P  +  r)  dP  +  «»*  ^  "  r)da 

A 


+  — - —  dx 
cos  v 


-Kj  -  cos 


(P  -  r)  do  -  2e  tan  9  cos  (P  -  r)  *P  *  eo*  ^  +  r)  40 


♦  dx 

cot  9 


-Kj  -  COS 


It.  - 


ja  -  (P  *  r) 

da  -  dvj. 

|a  -  (Pi  r) 

da  ♦  dv 

If  tht  d.t.mlnant  of  th.  •atrlx  of  co.fftd.ntt  In  Eg.  (3.3.*)  >  - 
not  t.ro.  th.  dariwati v.s  of  th.  «nhn<»nt  oay  b.  d.t.r.ln.d  onl,».ly  *'ong 
th.  curv.  y  -  y(x)  .  M0-...r,  If  th.  slop,  t.w  ca'vt  y  -  y<*)  »*  »««•» 

that  this  d.t.rwinant  do.*  vanish,  than  th.  d.rU.tlv.t  of  th.  onhiHX^ 
functions  along  th.  curv.  ar.  not  onigu.ly  d.t.fln*.  th.  curv.  I,  clUd 
a  ••charact.rittlc  " 

Th.  character  1st  let  -ay  b.  obtain*  by  tatting  th.  d.t.r«iMnt  0 
of  th.  co.ffid.ntt  In  Eg-  (3  3  *)  *«  *»r0 
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I 


0 


0  » 


This  determinant  is  equal  to: 


(3.3.8) 


D  - 


(3.3.9) 


Equating  the  first  of  these  determinants  to  zero,  leads  to  the  first  system 
of  characteristics. 


This  gives 


•) 

*2 


0  or  4o  tan  q  .  e(  a2  ■  0 


and 


4*  •  tan  (0  ♦  y) 


(3.3.10) 


•  tan  (0  -  y)  (3.3.11) 

where  y  is  the  angle  between  the  slip  planes  and  the  direction  of  greatest 
principal  stress,  tee  figs.  (2.2)  and  (2.3). 

These  are  the  characteristics  of  the  stress  field,  since  they  are 
a  consequence  of  the  stress  equation  fit*.  Since  a  system  of  two  r«el 
characteristic  directions  it  obtained,  the  original  equations  have  been  shown 
to  be  hyperbolic.  The  stress  field  characteristics  are  inclined  to  the 
■•axis  at  angles  (0  *  r) .  that  is,  at  the  tame  angles  at  the  lines  of  slip- 
Hence,  the  lines  of  slip  will  coincide  with  the  stress  field  characteristics 
on  the  ayplene. 


2S 


Equating  the  second  of  the  determinants  in  Lq.  (3.3.9)  to  zero  leads 
to  the  second  system  of  characteristics: 

C3 

c4 
or 

^  -  tan  (0  ±  r)  (3.3.12) 

^  -  -cot  (0  1  T)  (3.3.13) 

These  are  the  characteristics  of  the  velocity  field.  The  angle 
(0  +  r)  or  (0  -  r)  '»  used  according  to  whether  the  active  slip  line  makes 
an  angle  of  (0  *  r)  or  (0  -  y)  to  the  x-exis.  Thus,  one  of  the  velocity  field 
characteristics  coincides  with  the  active  slip  line  and  the  other  is 
orthogonal  to  it  as  is  shown  in  Fig.  3.1. 

The  characteristics,  or  lines  along  which  the  derivatives  of  the 
unknown  functions  may  be  discontinuous,  have  been  obtained  by  setting  the 
detcrminent  of  the  coefficients  in  lq.  (3-3.6)  equel  to  tero.  However,  for 
any  solution  to  exist  et  ell  for  the  derivetives,  the  equations  must  not  i>e 
inconsistent,  thet  is,  the  detcrminent  obtained  by  repleCing  any  column  of 
the  metric  of  coefficients  on  the  I  »f  t- hand  side  by  the  non  homogeneous  pert 
of  the  equations  must  also  vanish. 

Replacing  a  column  of  coefficients  on  the  left-hand  side  by  the 
column  on  the  right-hand  tide  of  lq.  (3.3.6)  and  setting  this  determinant 
equal  to  ser o,  yields  the  reietion  which  must  be  satisfied  along  the  charac- 
teristics.  As  before  the  reletiont  decouple  into  stress  field  end  velocity 
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field  enuetions. 


Replacing  the  third  column  on  the  left  of  Eq-  (3.3.6)  by  the  column 
on  the  right  and  setting  the  determinant  equal  to  zero  leads  to: 


K3  d3 


for  the  velocity  field. 

Evaluating  this  detcimlnant  yields: 


along  gj  -  tar.  (P  t  y)  i 


v  tan 


-  (p  *  r)J 


da  -  dv  »  0 


along  »  -cot  (0  t  y) 

01 


cot  j^a  -  (0  t  y)  j  da  ♦  dv  ■  0 


(3.3.14) 


(3.3.15) 


(3.3. .6) 


The  characteristic  equations  of  the  stress  field  art  obtained 
from  the  following  determinant: 


‘I 


‘I 


*3  CJ  d3 


»  0 


(3.3.17) 


tince  the  function  a  has  been  determined  from  the  velocity  field,  the 
coefficients  c,  and  Cj  times  the  derivatives  of  a  may  be  treated  ss  kn»*n 
quantities- 


It,  -  c, 


*2  ‘  C2 


-  0 


(3.3.  It) 
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This  yields: 


along 


along 


j*  -  tan  (p  +  r>;  K|  '  c,  g  -  0 

(3,3.19) 

2  -  tan  (P  -  r)  ;  *2  -  c2  t“  “  0 


Using  Eqs.  (2,3.7)  and  simplifying: 


cos 


dm 

IF 


do  ♦  2o  tan  ,  dP  +  7;~  ~4|  da 


(3.3.20) 


and 


cos 


d» 

IF 


*0  y 


COS  9  cos 


IF^D 


c 


(3.3.21) 


where  ($2)  and  (^)  are  evaluated  by  substituting  either  ^  *  tan  (P  ♦  r) 

or  2  •  tan  (P  -  r)  »"to  the  last  two  equations  of  the  system  (3. 3. ft). 

The  basic  equations  of  the  soil  medium  with  the  first  kinematic 
assumption  have  been  transformed  into  a  set  of  characteristics  and  a  set  of 
ordinary  differential  equations  to  be  integrated  along  each  of  the  character¬ 
istics. 


3.4  Character  i  st  i cs  of  the  fUcrp-Uructurel  Tiecy 

According  to  the  second  kinematic  relation,  the  velocity  vector  is 

di  reeled  along  the  active  slip  line  (see  Cq-  2.5-2).  from  Iqs  (3-2.7)  and 

(3.2  •)  the  velocity  components  v  and  v  are  given  es  functions  of  P  end  v: 

*  V 


v^  •  v  cos  (P  t  r) 

v^  •  V  Sin  (P  t  y) 


(3.4.1) 
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Using  these  relations  in  Eq.  (3.1.6)  satisfies  the  kinematic 
relation  of  the  macro- structural  theory  identically.  As  pointed  out  earlier, 
the  transformation  of  variables  (3.2.4),  (3.2.5),  and  (3.2.6)  will  satisfy 
Eq.  (3.1.3),  the  failure  criterion,  identically  as  well.  Substituting 
Eq.  (3.4.1)  into  the  equation  of  continuity,  (3.1.4)  provides  the  equation 
relating  v  and  6.  Using  this,  along  with  Eqs.  (3.1.1)  and  (3.1.2)  transformed 
to  the  new  variables  a  and  6,  one  obtains  the  basic  equations  of  the  macro- 
structural  theory. 


cos  (P  ♦  r)  ♦  ^  sin  (P  ♦  r)  ♦ 


2o  tan  9  p-  j  COS  (P  ♦  r)  *  t‘  COS  (P  -  y)  j  ♦  ^  j^sin  (P  ♦  y)  ♦  l*  sin  (P  -  y)  j| 


♦ 


A 

cot  « 


•  0 


(3  4.2) 


^  cos  (P  -  r)  *  sin  (P  -  r)  - 


2o  tan  9 


{*[ 


cot  (0  -  r)  4  t  tot  (G  ♦  r) 


]♦*[ 


(P  -  r)  *  I*  s>n  (0  r) 


J/ 


COS  9 


0 


(3  4.3) 


[* 


tin  (P  1  y) 


cot  (P  *  r) 


H 


s; 


cos  (P  t  r) 


*% 


tin  (9  t  y)  I  »  0  (3.4.4) 


where 


t 


v 

•j\o  tin  9 
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The  characteristics  of  these  equations  are  found  in  the  sane  way 


as  in  the 


xy-p lane , 


x  in  £qs. 


where 


continuum  theory. 

Consider  the  differential  rotations  along  a  curve  y  »  y(x)  in  the 


do  "  S  dx  +  f!  dy 

d0  •  dx  ♦  dy  (3.4,5) 

.  3v  dy  , 

dy  ’  S  dx  *  Ty  dy 


Using  these  relations  to  eliminate  the  derivatives  with  respect  to 
(3.4.2),  (3.4.3)  and  (3.4.4),  one  obtains  the  fol  lowing  equations: 


—  — 

-  - 

•1  b.  0 

£ 

*1 

*2  b2  0 

* 

m 

‘2 

!  b3  C3_ 

* 

*3 

(3-4.6) 


fc,  * 


b,  »  v 


dx  sin  (0  ♦  y)  -  dy  cos  (0  ♦  y) 

dx  sin  (0  -  y)  -  dy  cos  (0  -  y) 

2o  l*n  9  (e(  ♦  a?) 

*2o  tan  9  (a2  *  l*  a() 

j" dn  cos  (0  l  y)  ♦  dy  sin  (0  t  y)J 

dx  sin  (0  •  r)  *  dy  cos  (0  t  y) 


0^  ■  dx  sin  (0  I  y)  -  dy  cos  (0  t  y)  (3.4.7) 

Rj  *  cos  (0  ♦  y)  do  ♦  2o  tan  9  j  cos  (0  ♦  y)  ♦  (*  cos  10  -  y)  j  40  ♦ 

♦  *  <» 
cos  9 


3d 


-K^  ■  cos  (P  -  y)  da  -  2 a  tan  ©  [" cos  (P  -  y)  +  I*  cos  (fi  +  y)"j  dp  + 


+  8  dx 

cos  $ 

Kg  ■  v  sin  (P  +  y)  dP  +  cos  (P  ±  y)  dv 
From  Eq  (3.4.6), 


I2  * 

dy 


*! 

bl 

Kl 

K2 

b2 

&  . 

*2 

K2 

al 

bl 

?»y 

*1 

bl 

*2 

b2 

*2 

b2 

i 

(3.4.8) 


and 


-  3  3  dy 

^  =3 


(3.4.9) 


To  obtain  the  characteristics  of  the  stress  field,  the  determinant 
in  the  denominator  of  Eqs.  (3.4.8)  is  set  equal  to  aero.  This  gives: 

2c  tan  tp  J~ ?  a,  e2  +  t2  (e(2  +  ej2)j  •  0 
Using  Eqs.  (3.4.7),  this  simplifies  to: 


Jy  -(>•*•  t  sin  a)  sin  2P  t_cos  u  .  I  -  l 
***  (I  +  l2  sin  ©)  i.os  rP  *•  (J  4  sin  9) 


13.4.10) 


This  equation  gives  the  slope  of  the  character! stic  curve.  If  X  l»  the  angle 
between  the  x-axis  and  the  characteristic  cu»ve  at  any  point,  then 


*•"  X  *  £ 
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Let  r  be  the  angle  between  the  maximum  principal  stress  direction  and  the 
tangent  to  the  characteristic  curve-  Then, 

tan  f  «  tan  (X  -  0)  -  (3.4.11) 

Substituting  Eq  (3.4.10)  into  (3.4.11)  yields: 


or 


(3.4.12) 


(3.4.13) 


The  relation  between  the  slip  line  direction  and  the  characteristics  is 
shown  in  Fig.  3-2.  Equation  (3  4.10)  may  also  be  written  in  the  for* 


j£  -  t«n  (0  i  r)  (3.4.14) 

It  is  of  interest  to  consider  the  effect  of  the  factor  I  ■  —  ■* .  in 

<fla  sin  9 

Eqs  (3  4. 2)  to  (3.4.4)  for  various  values  of  (•  If  {  «  0,  this  implies 
that  v  ■  0.  This  is  the  static  case  (or  which 


tan  f  •  *  tan  y,  or  f  •  1  y 

and  the  stress  field  coincides  with  the  one  from  the  er**.inuu»  theory. 
If,  In  a  dynamic  cate  f  »  I , 


and 


v  •  \flo  si  n  9 


tan  a  ■  0,  or  f  •  0 
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for  this  casi  the  characteristic  curves  coincide  with  the  direction  of  maximum 


principal  stress.  When  the  velocity  exceeds  the  critical  speed  Jia  tin  9 
the  equations  become  elliptic,  and  the  method  of  characteristics  does  not 
apply.  The  existence  of  a  critical  speed  is  an  unexpected  consequence  of  the 
thcor,,  since  hyperbolic  equations  are  expected  in  dynamic  problems.  However, 
for  the  quasi-static  case,  the  difficulty  of  a  critical  speed  does  not  arise. 

The  differential  relations  to  he  satisfied  along  the  characteristics 
of  the  stress  field  are  obtained  by  setting  the  following  determinant  equal 
to  zero: 


tm 


(3.4.15) 


txpancing  this  determinant  and  using  Eqs.  (3.4.7)  and  (3.4.14),  one  obtains 
the  characteristic  equations  of  the  sfess  field. 


Along  ■  tan  (P  +  ?) 

dx  ('w+  St  cjdo  ♦  tan  9 

cosTE'TTT - - 


*  s.  c-w  +  <e- 


cos  q> 


h  * 


(3.4  16) 


Along  ■  tan  (p  -  f) 


dx 

cos  (P  -  fl  " 


r  1 

(c^y  +  cy)d 0  -  2o  tan  qr 

c  •  5,  t  +  r  (c  *  $.  t  )  Ml 
L  -V  1  V  *  V  1  -v'j 

cos  ft 


(3.4.17) 


where 


(3.4.  IB) 
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*  a 


T 


The  characteristics  of  the  velocity  field  are  obtained  by  setting 
Cj  of  Eq.  (3.4.9)  equal  to  zero.  This  gives: 


■gj  -  tan  (0  i  y) 


(3.4.19) 


Here  the  +  or  -  sign  is  used  depending  on  which  slip  line  family 
from  the  stress  field  is  active. 

Equation  (3.4.19)  shows  that  the  characteristics  of  the  velocity 
field  coincide  with  the  families  of  slip  lines.  This  is  to  be  eapectad  from 
the  basic  assumption  of  the  macro- structural  theory  which  required  that  the 
velocity  vector  at  any  point  coincides  with  the  direction  of  the  slip  line 
at  that  point. 

The  differential  equations  to  be  satisfied  along  the  velocity 
characteristics  are  obtained  by  setting  the  numerator  of  Eq.  (3.4.9)  to  zero. 


S  '  b3  $ 


(3.4.20) 


Using  Eqi.  (3.4.7),  this  becomes: 


Along  ^  •  ten  (0  ♦  r) 

<9  ♦  2e  sen  *  [  I  ♦  (c  ,  -  $u|  cot  »)]  #  ♦ 

ceVTjTVj "  rrr: 


"/cos  * 


(1.4.21) 


Along  •  ten  (0  -  y) 


dn 

cos  -  r) 


do  -  2«  ten  f  I  ♦  i  (c  .  +  S  2  cos  *)  d0  ♦  v  dv 

. . . . . -.-W .  jg . .i.iTii  i .  i .—*4 . . 

/cos  «  (3.4.22) 


± 


where 


c 


ul 


.  (P  -  r) 

cos  (g  +  r) 


ul 


sin  fp  +  r) 
cos  (p  +  r) 


c 


u2 


cos  (P  +  r) 
cos  (p  -  r) 


(3.4.23) 


3 . 5  Characteristics  in  the  Static  Theory 

it  will  be  useful  to  present  the  development  of  the  method  of 
characteristics  for  the  equations  governing  the  statics  of  soil  media*^ 
These  equations  are: 


.  3c  ck 
I  /  *  .  *y 


•« 


.  d*_  da 
V  *  l  ^  ’  ° 


2  2  2  2 
(o^  •  a  )  ♦  4  •  sin  ®  (%  ♦  Oy  +  2k  cot  *) 


Transforming  to  the  n«*  variables  c  and  P,  the  last  of  Cqt.  (3.5. 1)  Is 
satisfied  identically  and  the  following  two  equations  are  obtained: 

cos  (Pay)  ♦  ^  sin  (P  ♦  r)  ♦ 

♦  2o  tan  o  ^  tei  (»  +  r)  ♦  ^  *•"  (P  ♦  t)  j  ♦  •  ? 


(3S.1) 


(352) 


||  cot  (P  -  r)  +  ^  Sin  (p  -  r)  - 

-  2o  tan  «  j^j  cos  (P  -  r)  ♦  ^  sin  (P  -  r)  j  ♦  ~ "  •  0  '3S-3* 
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where 


A  -  X  sin  (0  -  r)  -  V  cos  (0  -  r) 

B  ■  -X  sin  (0  +  x)  +  X  cos  (0  +  r) 


Using  the  same  procedure  as  in  Section  3.3,  the  characteristics 
of  these  equations  are: 

-  tan  (0  +  r) ;  ^  (0  *  r)  (3.5.4) 

and  the  differential  equations  to  be  satisfied  along  the  characteristics  are: 

do  +  2o  tan  9  d0  ■  0 

(3.5.5) 

do  -  2o  tan  $  d8  •  0 

Equations  (3.5.5)  may  be  integrated  to  yield: 

C°j— ^  In  o  *  0  •  const  •  l(u) 

(3.5.*) 

In  0  ♦  0  •  const  •  *|(») 

Here,  a  function  of  o,  ainus  0  is  a  constant  on  Makers  of  one  faally  of 
characteristics,  and  a  function  of  o,  plus  0  It  a  constant  on  (Makers  of 
the  other  font ly.  Equations  (3.5.0)  nay  be  solved  for  0  and  the  function  of 


^y '*  l«  o  •  y  (|(u)  ♦  q(t)  ] 

0  •  y  l-l(u)  ♦  n(») ) 


(3.5.7) 


As  in  the  theory  of  plasticity,  special  cates  arise  when,  in  a  particular 

region  eeing  considered,  (  •**<  1  *r*  constants,  {  along  it  a  constant,  and 

.  .  .  15,10,17 

i)  alone  is  a  constant. 
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When  both  {  and  q  are  constant  functions  in  a  region,  Eq$.  (3.5.7) 
show  that  both  o  and  f  are  also  constants  throughout  the  region.  When  this 


is  the  case,  both  families  of  characteristics  are  straight  lines  and  the 
region  is  called  a  "region  of  constant  state"*^’*^. 

When  either  {  alone  or  t)  along  is  a  constant,  one  of  the 
characteristics  (3.5.4)  will  form  a  fan  of  straight  lines.  If  the  fan  is 
centered,  that  is,  if  the  straight  line  characteristics  intersect  at  a  point, 
the  region  is  called  a  "region  of  radial  shear"*. 


3 . 6  Transf  onset  ion  of  characteristics  to  Polar  Coordinates 

for  later  work  in  polar  coordinates  it  will  be  convenient  to  replace 
c:  and  0  b  -  similar  angles  Sr  and  b  measured  from  the  positive  direction  of  the 
radius  rector  r  ;nd  the  directions  of  the  velocity  vector  and  of  the  maximua 
principal  stress,  respectively.  The  angle  9  is  defined  as  the  counterclockwise 
angle  from  the  y>axis  to  the  radius  vector,  as  shown  In  fig.  (3.3) .  The 
fo!  'cuing  relations  may  be  easily  obtained  from  Fig.  (3.3): 


b  -  »  -  9  -  j 
5  •  (  -  8  ♦  y“ 


(3.b.  I) 


Po'ar  coordinates  will  be  used  later  when  In  a  aorticular  problem 
the  inactive*  family  of  slip  lines  consists  of  a  cente'ed  fan  of  straight  line 
character' st ics.  In  this  case  the  angle  b  in  polar  coordinates  has  a 
constant  value  equal  to: 


* 


That  is  to  say,  slip  actually  occurs  along  the  other  set  of  slip  lines- 
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tan  (p  +  r)  active 


(3.6.2) 


9  ■  +  r  for  ^  ■ 
dx 

or 

9  «  -  r  for  ^  »  tan  (p  -  r)  active 

In  polar  coordinates  the  velocity  field  characteristics  given  by 
(3.3.12)  and  (3.3.13)  become: 


d£ 

dr 


SSL. 8 

r 


de 

dr 


r 


(3.6.3) 


In  transforming  the  equations  to  be  satisfied  along  the  characteristics,  P, 
a,  da  and  dv  in  Eqs.  (3.3.  IS)  and  (3.3.16)  must  be  determined  fraa  the 
expressions  (3-6.1)  end  (3.6.2).  These  give: 


P  »  ♦  0  *  r 

°  *  *  +  *  *  T  (3.6.4) 

da  •  da  +  de 

dv  »  dv 

where  de  may  be  obtained  from  (3.6.3). 

Using  (3.6.4)  n  (3.3.  IS)  and  (3.3.16)  yields: 

tan  (5  *  2y!  ^40  *  -■"*  dr  j  •  ~  •  0 

(3.6. S) 

cot  (5  +  2y)  |" dO  ♦  ~ &  dr  j  ♦  ^  •  0 

where  '.he  •  or  4  sign  is  taken  depending  whether  the  active  slip  lines  arc 
inclines  at  an  angle  (P  ♦  y)  or  (P  -  r)  to  the  a-aals. 
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The  characteristics  in  polar  coordinates  for  the  second  kinematic 
assumption  may  be  easily  obtained  by  transforming  the  basic  velocity  Eq.  (3.4.4) 
to  polar  coordinates  and  then  finding  the  characteristics  f rom  the  transformed 
equation. 

Transforming  (3.4.4)  to  the  polar  coordinates  of  Fig.  3.2  yields: 


sin  ^  ±  ^  ^  t;n  <p  >  0 


(3.6.6) 


Again,  the  sign  in  the  second  term  is  selected  according  to  whether  the 
active  slip  lines  are  inclined  at  (0  +  r)  or  (6  -  y)  to  the  x-axis. 

Since  characteristics  are  lines  along  which  ^  or  may  be 
discontinuous,  and  along  which  the  relation 


dv  "  t?  dr  *  ^  * 


(3.6.7) 


holds,  the  system  (3.6.6)  and  (3.6.7)  may  be  solved  in  the  usual  way  to  yield 
the  differential  equations  of  the  character) sties. 


*  »  ,  ££L 

dr  r 


(3.6.8) 


Integrating  (3.6.8)  gives  the  equation  of  the  velocity  field 
characteristics  and  active  family  of  slip  lines. 


r  ■  r  e 


1  8  tan  9 


(3.6.8) 


The  differential  equation  to  be  satisfied  air  .g  the  characteristics 


(3.6.9)  Is: 


or  upon  integration: 


~  -  3  tan  *  dd 


v  .  v 
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(3.8.10) 


(3.6.11) 


4.  NUMERICAL  SOLUTION  OF  EQUATIONS 


4. I  Finite  Difference  Method  of  Solution 

In  general,  the  solution  of  the  differential  equations  to  be  satis¬ 
fied  along  the  cha-acter i st ics  can  not  be.  obtained  analytically.  However,  a 
solution  can  be  obtained  numerically  by  transforming  the  equations  to  charac¬ 
teristic  coordinates,  and  then  converting  the  differential  equations  Into 
difference  equations  along  the  characteristics. 

The  region  in  which  the  solution  is  to  be  obtained  Is  divided  into 
a  network  by  the  uz-  characteri Stic  coordinates  as  shown  in  Fig.  4.1.  The 
solution  proceeds  in  a  step-bystep  fashion  by  obtaining  the  values  of  the 
unknowns  at  the  point  &k  t  after  the  values  of  the  preceding  points, 6^  ^ 
and  t  have  been  obtained.  At  the  start  of  the  computations,  points 

S,  .  ,  end  b,  ,  ,  correspond  to  points  on  the  Initial  curve  for  which 
k,l-l  k- 1  . 1 

initial  values  are  prescribed. 

Let  u(x,y)  be  the  coordinate  measured  along  one  family  of  charac¬ 
teristics  and  a (x,y )  the  coordinate  along  the  other  family.  In  the  ux- 
plane 


and 


Solving  for  dx  and  dy 


du  •  Q  4*  ♦  dy 

(4.1.1) 

Hti  -*•(£)) 
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where 


du  dz  bu  dz 
lx  |y  ly  lx 


Since  x  -  x(u,z)  end  y  ■  y(u,z) 


dx  •  ^  du 


dz 


(4.1.3) 


“>'  "  fc  dU  + 

Equations  (4.1.2)  and  (4.1.3)  yiald  tna  fol lowing  transformation  relations. 

it -*lf  >  *--•* 

(*.!.«) 

if-'t* 

|f  for  any  function  f,  f  ■  f(u,z)  ,  then: 

s -[«•«] 

(4. I . &) 

perivatives  with  respect  to  the  characteristic  coordinates  u,z  are  epproxi  - 
mated  by  the  following  finite-difference  expressions: 
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,Yi ;  Vi,i 

Au 


d*  wXk,l  '  xk,i-l 

5z  Az 


ksyk  ,i :  vi.i 

du  &u 


if 


^yk,i  •  yi,M 
Az 


(4.1.6) 


M  -V*  I 

du  Au 


df  g  ^k,l 

5z  Az 


-  f 


k.i-l 


4.2  Equations  In  flnlt.~Dlfftr.nce  form 

Along  th.  ch.r.ct.r 1  it ic  u  «  u(x,y)  -  const., 


14.1.1) 


If  th.  v.locity  fiald  character I st Ic  glvan  by  (3.3.12)  Is  Masd  th.  on.  along 
which  u  -  u(x.y)  -  const.,  Eg.  (4.2.1)  baconas 


(M  r)  ♦  ^  ala  (?  4  r)m  0  (4.1.1) 

Similarly,  along  i  ■z(x.y)  ■  const.,  using  Eg.  (3.3.13) 

r)  - cm  (|i  s  r)  «  j  (4.i.» 

Using  (4.1.4)  In  Eg*.  (4.2.2)  and  (4.2.3)  ylatds 

do  -  j*,  (p  «  r)  .  |5  »ta  (A  i  f)  -  0 

(4.2.4) 

ds  *  |J  sla  ($  4  f)  ♦  SM  (fi  4  f)  -  0 


i 

t 

i 

i 

I 

♦ 

i 

i 


j 

i 

( 
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These  are  the  characteristics  of  the  velocity  field  In  the  uz -piano. 


In  a  like  manner,  using  (4.1.4),  (4.1.5)  and  (4.2.4)  in  Eqs.  (3.3.15) 
and  (3.3.16),  one  obtains  the  equations  to  be  satisfied  along  the  character¬ 
istics  in  the  form 


vten  [«-(*tr)]  8-fe  -  0 
v  cot  -  (0  *  Y)  2  $  ♦  jj  -0 

Equations  (4.2.4)  and  (4.2.S)  may  be  expressed  In  f Inlte-dlfference  form  by  making 
u  e  of  the  difference  quotients  (4.1.*). 


V.,- 

[*k.l  ’  %J-I 

]  •••  (  *k.4-»  *  r  /  *  0 

(4.1.4) 

k.i "  Vi,<]  ,,B  r  )  ♦ 

[“k.i  '  Vl,4. 

|  <w  )  •  « 

end 


tan 


cot 


=  -JuirJ 

b,*- 1 

4 


0 


(*.1.7) 


The  finite -difference  form  of  the  aqua t lens  of  the  velocity  field  character¬ 
istics,  fqe.  (3.4.3),  In  paler  seerdlnetae  Is 


(V«  -  *k 


Vl'  ek 


(4.2.8) 


The  aquations  along  tha  characteristics  taka  the  finite-difference  for« 

tan  fak  *  2/)  ^hiJ  .  “k,i-i  -  cot  -(VH,Jk"^id) 


sk-i.«  * l7>  ^k,«  -  ;k-i.<  ♦  *('l,‘lv^.l,<) 


S.  PRESSURE  LOADING  ON  BURIED  STRUCTURES 


5.1  Stress  Field  Soiut Ion 

The  theory  presented  in  the  previous  chepters  will  be  epptled  to 
the  problem  of  determining  the  intensity  of  pressure  ecJng  on  e  structure 
buried  in  e  grenuler  soli  when  the  ground  surface  is  subjected  to  an  air- 
induced  over-pressure  and  the  structure  is  ok  the  verge  of  collapse.  For 
this  case,  the  oressure  applied  to  the  structure  corresponds  to  e  Measure  of 
its  ultimate  strength. 

The  structure  considered  is  a  plate  having  one  side  much  longer 
than  tha  other  so  that  tne  two-dimensional  state  of  plane  strain  will  prevail 
at  regions  away  from  the  ends  of  the  long  side.  Taking  the  xy-plane  perpendicu¬ 
lar  to  the  long  axis  of  the  plate,  the  slip  planes  will  form  as  slip  lines  on 
the  xy-planc. 

The  intensity  of  pressure,  pQ,  applied  to  the  ground  surface  Is 
assumed  to  be  uniformly  distributed  over  a  large  area,  and  directed  vertically. 
Since  there  Is  no  shearing  stress  at  the  ground  surface,  the  vertical  stress 
in  the  soil  at  the  surface  Is  a  principal  stress  ague!  to  p  .  The  structure 

9 

>s  assumed  symmetrical,  so  by  virtue  of  ryemetry  the  vertical  stress  along 

the  axis  of  symmetry  AB  in  Fig.  5.1  it  also  a  principal  stress. 

Ir  general,  at  in  tha  Theory  of  Plasticity,  tha  stress  field  equations 

can  not  be  ipplled  directly  to  a  specific  problem  obtain  tha  required  failure 

praisure.  However,  a  solution  It  obtained  by  finding  a  proper  stres*  field 

for  which  there  exists  a  compatible  velocity  field  satisfying  alt  velocity 
ia 

boundary  conditions 

The  families  of  stress  characteristics  or  slip  linos  have  certain 
geometric  properties  which  fectUtata  tha  construction  of  s  strost  field 


from  given  conditions  at  the  boundary.  These  oroperties  were  first  presented 
!S  !9 

by  Hencky  and  Prandtl  for  the  Mises  Theory  of  unrestricted  plastic 
flow.  One  of  the  most  important  of  these  geometric  properties  is  that  the 
angle  formed  by  the  tangents  of  two  slip  lines  of  one  family  at  the  points 
of  intersection  with  a  slip  line  of  the  other  family  is  a  constant  no  matter 
what  member  of  the  other  family  is  taken. 

From  Eg.  (3.3.9),  it  is  seen  that  the  stress  and  velocity  fields 
■t'.j y  be  considered  separately.  Additional  properties  of  the  slip  lines  are 
given  by  the  relations  derived  for  the  static  theory  of  Section  3.5  when  the 
quasi-static  problem  is  considered  in  which  inertia  forces  are  not  considered. 
Since  the  intensity  of  applied  over-pressure  is  a  constant,  a  region 
of  constant  state  consisting  of  two  families  of  straight  slip  lines  forms  at 
the  ground  surface-  Similarly,  assuming  a  uniform  pressure  applied  to  the 
structure,  we  have  a  region  of  constant  state  adjacent  to  the  structure 

The  horiiontal  principal  stress  in  the  soil  at  the  ground  surface 
is  taken  to  be  the  maximum  principal  stress  as  a  result  of  arching  in  the  soil- 
Terxaghi  ^  has  verified  this  behavior  experimentally.  With  the  x-axis  vertical, 
the  value  o'  in  the  constant  state  region  at  the  surface  is  either  tx/2-  The 
d. section  of  m«ximur  principal  stress  at  the  structure  is  vertical  so  that 
the  angle  ft  undergoes  a  total  change  of  */2  and  is  equal  to  xero  in  the  constant 
st„.c  region  adjacent  to  the  structure 

figure  5.2  shows  the  network  sf  characteristics  for  the  stress  field 
solution  of  the  guati-stetic  probelm.  The  stress  field  is  di'  ide.d  into 
four  regions,  the  two  c  instant  state  regions,  COC*  and  COO*  and  the  ceniered 
fans  OCD  and  0C’0‘  .  In  syiematr ic.i  problems  both  families  of  slip  lines  will 
b''  active  11  ,  se  that  along  kht  limiting  active  slip  line  CO,  dy/dx  •  tan(P-r) 


and  along  C'D',  dy/dx  ■  tan(0+r),  The  u  and  *  families  of  characteristics  are 
given  in  region  00D‘  by  the  differential  equations  dy/dx  ■  ten(3±y).  In  this 
region  'he  angle  3  -  const.  ■  0,  since  the  direction  of  the  maximum  principal 
stress  coincides  with  the  direction  of  the  x-axIs.  Then  tire  equations  of 
the  slip  lines  ere  given  by 


>  »  -  x  tan  x  *  c| 
y  ■  x  tan  y  +  e. 


(u-famlly) 

(z-famlly) 


where  c,  and  c?  are  arbitrary  constants  and  the  u-famlly  is  active. 

The  characteristic  for  region  COC'  (3  -  const.  -  n/2 )  are  obtained 


in  a  s i ml iar  wav  to  be 


y  •  x  cot  y  +  Cj 
y  •  -  x  cot  y  +  c. 


(u-fanily) 

(z-family) 


(5.1.2) 


In  region  0C9  it  is  convenient  to  express  the  character istlc  net 
by  a  system  of  plane  polar  coordinates  with  9  measured  countarclocfcwttn  as 
shown  in  Fig.  5.2.  Mere,  the  direction  of  the  Maximum  principal  stress  Is 
given  by  the  angle  0  which  Is  measured  with  respect  to  the  radius  vector 
from  point  0  taking  counterclockwise  positive.  Since  the  redll  fro*  point  0 
ere  members  of  the  Inoctlve  family  of  si  ip  lines  In  region  0C9,  the  angle  0 
Is  eonstent  in  this  region  (cf.  Eqs.  >.6-2)  'n  pnlar  coordinates  the 
differential  equations  of  the  stress  field  chereeterlttici  in  region  0C9  ere 


^  (u-fasl  ly) 


f- 


(5.1.3) 


(■-family) 


They  yield 


r  .  rje-«Mr>tan<P  (u-fami  ,y) 


{5.1.4) 


e  -  e 


(z-femi ly) 


The  u-family  of  active  slip  lines  form  a  family  of  logarithmic  spirals  given 
ty  tne  first  of  Eqs.  (5.1.4). 

The  relations  to  be  satisfied  along  the  characteristics  given  in 
Section  3.5  are 


where 


ii  -  &  "  (  •  const,  on  u-faml ly 


a  +  8  »  n  *  const,  on  z-tamily 


(5.1.5) 


„  -  ,  .  o  +  k  COt  T>  i 

-  .  cot  9  in  » 

0  2  [l  +  sin  V  cos  20 


,%  the  constant  state  region  000’,  both  |  and  g  are  constants.  Thus.  this  region 
transforms  info  a  finals  point  on  the  I  Ine  l  •  g  on  the  Ig-charactarlstlr. 
plane  shcmn  in  Fig.  5-3. 

On  the  right-hand  side  of  the  axis  of  symmetry  AS,  the  maximum 
principal  stress  direction  changes  from  0  ta  90  degrees  with  respect  to  the 
x-nxis  by  a  countarclockwlsa  (positive)  rotation,  -ao  tne  angle  0  Ir.  region 
ACC  Is  ♦  e/2.  On  the  left-hend  side  of  A*,  the  w  dmum  print  I  net  St.*#** 
direction  changes  by  e  clockwise  (negative)  rotetlon  of  90  degrees,  so  the 
value  of  B  In  region  AOC'  is  -  */2.  in  general,  the  propei  *lg.>  of  the 
tngle  0  must  bh  used  In  determining  t*#  correct  valiaet  of  the  'tresa  com¬ 
ponents  . 


| 
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Regions  ADC  and  AOC1  are  shewn  In  Fig.  5.3  as  points  on  the  iq-plane. 


The  centered  fans  OCD  and  OC'D'  transform  Into  the  lines  {  ■  const,  end  T|  ■  const, 
connecting  the  constant  state  regions.  Although  the  region  COC1  Is  transformed 
into  two  points  on  the  characteristic  plane,  only  one  need  he  considered  since 
each  point  yields  the  same  values  of  stress  components. 

The  applied  pressure  po  required  to  develop  the  assueed  stress 
fie'd  may  now  be  determined.  Using  the  first  of  Eqs.  (S.I.S)  to  determine 
the  value  of  t  in  regions  AOC  and  MO  In  terms  of  the  pressures  pQ  end  q, 
it  may  be  seen  from  Fig.  5.2  and  Fig.  5.3  that  these  values  must  be  equal . 

By  equating  the  values  of  and  the  pressure  pQ  Is  obtained  as  a 
function  of  the  ultimate  plate  resistance  q. 

In  region  AOC,  o  •  p  and  0  ■  */2.  From  the  first  of  Iqs.  (S.I.S) 

K  9 


lAOC  *  'am 


In  region  IOC,  p  •  q  and  0  •  0.  Then 


<Mt  *  <MC 


“F  ’fb^]  •  I  •  T IHW] 


"‘fi'rrai '  Fife?] 


*  «  tan  f 


(S.1.7) 


(>.!.*) 


Then 


-0  •  <- +  - 


x  tan  7 


-  k  cot  9 


(5.1.10) 


For  9  «  30*.  k  ■  0  (Ideal  granular  material) 


P0  -  q| 


fl  -  tin  9~)e: 
[l  +  tin  9J 


x  tan  9 


or 


|>o  •  2.03q 

Expresting  the  soil  cohesion  k  by  the  relation 


Ck< 


Eq.  (5.1 .10)  becomes 


P0  -  «[<» 


♦  Ckcot  9)f(9)  -  C^cot  t 


(5.1.11) 


•(here 


C 


k 


.k 

q 


f(9) 


s  tan 


* 


It  should  be  noted  that  this  solution  correspond*  to  a  particular 
depth  of  burial  depending  on  the  span  length  l  end  the  angle  of  Internal 
friction  of  the  soil  9.  evaluating  the  arbitrary  constant  r^  In  the  first  of 
tqs.  (5.1.4),  the  equation  of  ths  line  CO  is  obtained 


r  •  00  a 


-(0  -  r  ♦  fitan  9 


(5.1.12) 


10 


The  length  OC  is  found  by  setting  6  *  +  T 


OC  -  00  e  2 


-  *  tan  <P 


(5.1.13) 


From  this  the  length  CC‘  and  the  depth  H  may  be  found  in  terms  of  the  span  L. 


CC  -  L  cot  r  e  2 


-  t  tan  <P 


(5.1.14) 


H  -  L 


1  -  ^  tan  9  .  I 

2  +  j  cotr 


(5.1.15) 


for  9  -  30* 

r 

Then 

CC‘  -  O.TOL 
and 

H  •  107l 

it  should  be  remembered  that  the  solution  obtained  hare  It  valid  for  only  on* 
depth  if  the  span  of  the  structure  and  9  It  given. 

5.2  Solution  for  Arbitrary  freftthi 

The  tolutlon  for  the  stress  field  obtained  In  lection  5.1  It  based 
on  the  ettunption  that  the  entire  region  bounded  by  the  limiting  ti  p  Unas 
CO  and  CO'  It  In  a  ttete  of  plattlc  equilibrium,  lltlng  thli  atsiatptlon  and 
the  ttrett  field  of  Fig.  $.2  leads  to  a  particular  ratio  of  depth  to  span 
length  for  given  toil  properties.  The  velocity  field  correepondlng  to  tMt 
street  field  will  be  given  leter. 


St 


When  the  depth  of  burial  is  greater  or  less  than  the  "critical"  depth 
yiven  by  £q.  (5.1.15),  the  failure  region  bounded  by  the  limiting  slip  lines 
may  contain  zones  which  remain  essentially  in  the  elastic  state.  In  these 
elastic  regions  the  distribution  of  stresses  is  unknown  since  the  basic 
equations  of  the  analysis  do  not  apply.  However,  if  an  elastic  zone  forms  on 
a  boundary  of  the  failure  region,  it  may  be  possible  to  determine  the  total 
pressure  transmitted  by  the  elastic  region  to  the  interior  and  to  obtain  a 
solution  without  knowing  the  stress  distribution  completely.  Such  an  approach 
was  used  by  Hill,  Lee  and  Tupper  for  the  problem  of  a  ductile  material 

compressed  between  rigid  plates. 

Solutions  for  depths  less  than  the  "critical"  may  be  obtained  by 

assuming  that  regions  COC  and  00D1  are  not  fuily  plastic  and  that  the  angle 

0  in  the  plastic  regions  OCO  and  0C“B;  undergoes  a  change  of  lass  than  90 

degrees  m  going  from  00  to  0C.  If  and  *00  are  the  values  of  P  along  0C 

and  00  and  p  and  q  are  taken  to  be  equal  to  a  on  0C  and  00,  then  the  failure 
o  * 

pressure  p  Is  given  in  terms  of  q  by 
o 


{«■ 


♦  C^cot  »} 


1  ♦  sin  9  cos  2P, 


I  ♦  sin  9  cos 


5$  [•’  “  ■ V]  •  V-  •} 


(S’  ') 


•Aer* 


k 

9 


The  depth  of  burial  H  will  now  be  found  In  tanas  of  P^  and  P^. 


The  equation  of  the  curve  CO  Is  given  by 


4 


rco  -  OD  e*^0  +  I  ‘  r  '  Poo^tan  V 


^inLe  e0C  r  ~  4  :^nr '  Eq ■  (5-2.2)  becomes 


Oli  -  OD  e'^OC  ‘  ^OD^tan  ^ 


(5.2.2 


(5.2.3) 


OD  «  ^  •  —  ■  1 - 

2  sin (r  +  ft 


OA  -  OC  sin(r  +  90(.  -  j) 


the  depth  H  and  the  length  CC'  are  given  by 


"  ■  5  ;">*<'  *  »o„)  *  ;>  '  B°°h‘r  ’]  (5.2-4) 


■  ■■  jTTVtr  I  pj  •  V'-  *' 


(5.2.5) 


•*(n  ^QC  '  "/J  and  ’'no  *  °*  Eq*  (5  2-4)  •«<»  (5.2.5)  reduce 


•  ?  [cot  r+  e*  2  tan*] 


CC*  «  L  cot  r  e’  2  ta"  * 


(5.2.6) 


(5.2.7) 


These  are  the  same  as  Eqs.  (5.1.14)  and  (5.1. IS)  given  in  Section  5.1. 


.H 
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Equations  (5.2.!)  and  (5.2.4)  yield  a  range  of  failure  pressures 


I 

I 


for  any  particular  depth  depending  on  the  values  of  Bq^  and  Bqq  used.  The 
proper  values  of  B^  and  Bqq  are  the  ones  which  yield  the  smallest  failure 
pressure  at  a  given  depth.  Figure  5.4  shows  the  configuration  of  the  failure 
region  for  the  critical  depth  and  for  three  other  cases  in  which  the  depth  is 
less  than  critical.  These  cases  are  for  a  cohesionless  material  with  an 
angle  of  internal  friction  equal  to  30  degrees. 

When  the  depth  of  burial  is  much  greater  than  the  critical  depth, 
the  state  of  stress  at  the  ground  surface  will  be  affected  only  slightly  by 
the  arching  of  the  soil  above  the  yielding  structure.  It  has  been  shown 
exper imental I y  7  that  the  pressure  of  a  yielding  structure  in  a  granular 
medium  has  little  or  no  effect  on  the  state  of  stress  at  the  ground  surface 
when  the  depth  of  burial  is  on  the  order  of  2  to  3  times  the  width  of  the 
yielding  element.  If  the  state  of  stress  at  the  surface  !s  not  a'fected  by 
tne  yielding  structure,  then  it  may  be  concluded  that  failure  as  defined 
herein  will  not  occur.  The  structure  will  be  safe  provided  that  it  has 
sufficient  deformation  capacity  to  absorb  the  energy  of  the  initial  impulse 
applied  by  the  pressure  wave. 

Solutions  for  depths  greeter  than  the  critical  mav  be  obtained  by 
considering  the  stress  field  at  the  critical  depth  to  be  extended  by  a  region 
consisting  of  elastic  tones  and  a  centered  fan  In  sdiich  the  angle  8  undergoes 
a  change  of  less  than  90*.  As  the  depth  increases,  the  change  m  8  approaches 
90*  and  B  at  the  ground  surface  approaches  0“ ,  the  value  corresponding  to  the 
state  of  stress  existing  in  the  soil  if  no  yielding  structure  was  present. 

A  stress  field  for  a  depth  greater  than  the  critical  is  shot**  in 
Fig.  5.s(a).  Considering  the  regions  to  the  right  of  tho  axis  of  symmetry, 
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the  failure  region  consists  of  the  cr i t ical -depth  field,  regions  SAC  and  566, 
and  the  centered  fan  OCD  in  which  g  =  const. 

Regions  OAC  and  OBD  are  taken  to  be  in  the  elastic  state,  that  is, 
they  move  as  rigid  bodies.  As  in  the  case  of  the  solutions  at  depths  less 
than  the  critical,  the  velocity  field  associated  with  this  stress  field  will 
develop  if  the  structure  moves  vertically  downward  with  a  uniformly  distributed 
velocity  so  that  the  clastic  regions  may  move  downward  without  deforming.  The 
'uilure  pressure  obtained  on  the  basis  of  this  assumption  should  be  a  minimum 
since  any  other  motion  of  the  structure  requires  the  elastic  regions  to  become 
partially  or  completely  plastic  and,  thus,  more  work  must  be  done  by  the 
external  loads  in  defo.ming  them. 

<n  the  centered  fan  6c6,  the  angle  9  changes  by  an  amount  varying 
from  0"  to  90*  depending  on  the  depth.  When  the  changa  in  9  is  less  then 
90°,  the  point  on  the  (n*plsne  representing  region  666  falls  between  the 
9  •  s/2  and  9  •  0  lints  as  shown  In  Fig  5.5(b) 

Figure  5.6  shows  the  stress  * i « I d  for  the  critical  depth  and  thrae 
other  cases  at  a  greater  depth  for  a  coheslonless  material  with  T  •  JO*.  It 
may  bt  observed  from  the  values  of  fellure  pressure  given  that  above  the  critical 
depth  the  pressure  required  for  failure  increases  very  rapidly.  This  behavior 
should  be  expected  on  the  basis  of  tha  experimental  results  mentioned  previously 
The  effect  of  angle  9  on  feiture  pressure  at  critical  derth  H  Is 
shown  «n  Fig.  5  7.  The  sensitivity  of  failure  <  rassure  to  varlatloM  of  9  is 
apparent  here  A  plot  of  failure  presture  versus  depth  is  given  in  fig.  5  6 
for  "5  ■  30*  and  k  ■  0-  Here  the  targe  increase  In  failure  pressure  for  depth: 
greet f  then  the  critical  is  clearly  shown. 


5  3  Veiacitv  Field  Solution  for  the  Continuum  Theory 

The  equations  of  the  velocity  characteristics  and  the  differential 
relations  to  be  satisfied  along  them  for  the  quasi-static  incompressible  case 
are  (cf.  Section  3.3): 

u  ■  const •  fami ly 

P-  ■  tan(0  t  r ) ;  v  tan  [a  -  (8  ±  r)}  da  -  dv  ■  0  (5.3.1) 

dx 

z  -  const,  fami ly 

•  .  cot (8  t  r);  v  cot  (a  -  (6  i  y)}  da  +  dv  •  0  (5.3.2) 

dx 

for  the  regions  to  the  right  of  the  axis  of  syenetry  AS,  the  active 
family  of  slip  lines  is  given  by  the  equation 

^  •  t»fl(?  t  T)  (5-3.3) 

OX 

Thus,  one  fami ly  of  velocity  characteristics  is  coincident  with  the 
active  family  of  slip  lines  and  the  other  family  Is  orthogonal. 

The  velocity  field  characteristics  are  shown  In  Fig.  S.»  in  the  xy- 
plane  and  In  Fig.  5.10  In  the  ui-charecterlst le  plane.  In  region  010  the 
angle  8  maintains  the  constant  velue  determined  by  the  stress  solution  to  >e 

8-0 

Then,  In  this  region,  along  thr  lines  u  -  const. 

Ij-iany;  t#n(a  ♦  y)<a  -  ^  -  0  (5.1.4) 
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and  along  z  •  const . 

^  -  cot  r;  cot(a  +  r)dar  +  ^  *  0  (5.3.5) 

Integrating  the  equations  to  be  satisfied  along  the  characteristics  yields 

along  u  *  const.;  v  cos  (a  +r)  *  Cu(u)  (5.3.6) 

along  z  *  const . ;  v  sin  (a  ♦  r )  ■  C^(z)  (5.3.7) 

The  initial  velocity  distribution  will  be  specified  at  the  boundary  formed  by 
the  yielding  structure. 

It  is  assumad  that  the  soil  particles  immediately  adjacent  to  the 
yieiding  structure  will  move  in  a  vertical  direction,  that  It,  a  •  0.  Also, 
by  virtue  of  symmetry,  the  soil  particles  must  move  vertically  (a  •  0)  along 
the  axis  of  symmetry  AS.  These  boundary  conditions  are  given  by 

on  M  v  •  v  (y)  (5.3.6) 

o 

a  •  o 

on  At  O  •  3  (5-3.6) 

Since  both  a  and  v  are  known  along  the  non-cherecterlstic  68,  the 
solution  me,  be  obtained  In  the  triangular  region  formed  by  60  and  the 
characteristic'  66j  and  DBj  through  the  points  6  and  6.  6e*t,  tlx.  solution 
may  be  obtained  ir  the  region  66^1  since  3  l»  along  the  non-Character ist ic 

6 C. 

At  this  stag*.  Cv*>dltion*  in  the  curvilinear  region  0C6  mutt  be 
conside~«d- 
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In  region  0C0  it  is  convenient  to  express  the  equations  in  polar 
coordinates  r,6  with  r  measured  from  the  origin  of  the  xy-coordlnates  and  9 
measured  from  the  positive  y-axis,  positive  counterclockwise.  In  the  new 
coordinates,  0  and  Of  represent  the  angles  between  the  positive  direction 
of  the  radius  vector  ard  the  directions  of  the  maximum  principal  stress  and 
the  velocity  vector  Vj  respectively,  taking  counterclockwise  positive.  By 
using  oelar  coordinates  in  region  QCD,  the  angle  3  has  a  constant  value 
equal  to 

3  -  x  -  T  (5.3.10) 


In  Section  3.6  it  was  shown  that  the  angles  3  and  a  In  polar  coordinates 
are  related  to  0  and  a  by 


5  -  0  -  9  •  | 
a  •  a  -  8  ♦  y* 


(5.3.11) 


The  characteristics  In  polar  coordinates  and  the  relations  to  be  satisfied 
along  them  for  dy/dx  ■  tan (0  -y)  active  are  (cf.  tqs-  3.6.3  and  3.6.5} 
along  u  *  const . 


A  .  s-ttLf 

4r  r 


(S.J  II) 


.-(5 .  ir)[«  -  asilt]  -  *  -  o 


; j.3.i3) 


along  x  »  const - 


A.smlI 

4r  r 


(6.3. }«) 


5* 


cot(3  +  2r>!  da  +  +  *  .  o 

I  r  J  * 


(5.3.15) 


The  Eqs.  (5.3.12)  and  (5.3.14)  may  he  ntegrated  to  yield  the 
characteristics 


u  *  const . 


z  •  const . 


r  -  r  e"{0  +  2  "  r)tan  9 
o 


r.r*(9+i  -rJwfP 
o 


where  r^  is  measured  along  00 . 

Along  the  limiting  slip  line  CD,  the  normal  component  of  velocity 
vanishes'1  so  that  the  velocity  vector  at  any  point  on  CD  has  a  direction 
given  by  the  tangent  to  CD  at  that  point. 

Then  along  CO  the  angle  a  is 


ra  »  2s  - 


(5. 3. 16) 


Substituting  this  value  of  a  into  (5-3.13),  the  aquat  on  which 
holds  along  the  u  *  const  line  CD,  yields 

JZ  CO.t_f_dr1  dv  n 


tin (2k) 


dr*  - 


r  j  v 

Since  ton (2n)  *  0,  this  equation  shows  that  along  CO 


or 


dvco  -  0 


WCD  "  coo,t' 


The  value  of  a  in  xy-coordinates  at  point  D  It 


(5.3.17) 


°P  *->• 
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Substituting  this  into  (5.3.5),  the  differentia'  equation  to  be 
satisfied  alonu  the  z  =  const,  line  DE ,  yields 


^DE  =  0 


(5.3.18) 


\y  % 


,  t 


.tDE  =  const .  »  -  r 


Equation  (5.3.18)  «hows  that  the  angle  a  undergoes  an  abrupt  change 
,  ’ong  the  characteristic  OE  Substituting  the  value  of  O.  given  by  (5.3.13) 
into  (5  3.7)  shows  that  v^E  is  arbitrary.  The  value  of  v^E  to  be  used  along 
OE  in  extending  the  solution  is  obtained  from  the  condition  that  the  components 
of  the  velocity  vector  normal  to  OE  must  be  continuous  v;hen  the  direction  of 
the  velocity  undergoes  an  abrupt  change  across  OE . 

The  value  of  v  along  CO  may  now  be  obtained  from  the  initial  velocity 

at  point  0 


vco  *  cos  r 


(5.3.19) 


<r 


Both  '1  and  v  are  now  known  along  the  character ist ics  CO  and  OE 

and  the  solution  may  be  extended  into  the  region  OE. E bounded  by  the  Initial 

characterise ics  OE  and  DC  and  the  characterise ics  Ef^  and  CEj-  The  solution 
it,  nc*:  obtained  in  -eoion  EEjf  in  me  same  way  as  region  M(£.  This  method 

is  continued  until  the  solution  is  obtained  in  the  region  BOCHj- 

A  method  of  finite-differences  Is  used  in  obtaining  the  solution  In 
the  curvilinear  region  OCD  since  the  equations  to  be  satisfied  along  the 
characteristics  cat.  not  be  integrated  a s  they  could  it.  region  0*0.  for  the 
f inile-differtnee  solution  the  region  tt  divided  Into  i  network  of  u  and  z 
ch«r»c;erei'.t  ics  and  the  node  points  are  numbered  by  the  coordinates  k  and  t 
taking  i  along  the  u  *  conit  lines  and  k  along  the  t  •  const.  I  mas  • 


* 


O 


Expressed  in  finite-difference  form,  the  equations  along  the 
characteristic  become 

along  u  const 


la"(5k.i->  +  2r){ak,i  '  \.t- 


-  cot  <P 


along  z  “  const. 


[rk.4  *  rk.4-l 

1 

rvk.e  ■  vk.i-il 

rk,4-l 

/ 

L  vk,4-l  -1 

(5.3 

rrk ./  •  vi./i i . 

[Vi  ■  Vi.i] 

L  Vi. 4  - 

j T 

L  vk-1.4  -1 

«=  0 


(5.3.20 


Knowing  the  values  ef  a  and  v  at  tnt  points  k,4-l  and  k-1,4  tha 
above  algebraic  equations  My  da  solved  to  yield  a  and  v  at  tha  point  k,4 . 

The  solution  for  tha  values  of  <3  and  v  at  node  points  on  tha 
line  00  are  obtained  by  solving  (5.3 .6)  and  (5-3.21)  simultaneously.  Using 
(5,3. 1 1) ,  tha  equations  to  be  solved  for  points  on  00  are 


„  .  Vi-, i, 

k.i  cos(ok  |  ♦  r 


♦  xt 


(5,1.22) 


and 


t.4  "Vi.i  W<3k-1 .1  *  *>{Sk.4  *  Vi.i  *  '*  }  4  Vi. 4 


(5.3.33) 


where 


Vl.l  k-1,4 


k-  I  4 


1 

2 
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The  solution  for  OAC  is  carried  out  in  essentially  the  same  way. 


5  4  Velocity  Field  for  the  H/.tro-sf ructvral  Theory 

The  velocity  field  solution  obtained  with  the  second  kinematic 
relation  will  bo  presented  for  the  corresponding  stress  field  solutions  or 
Sections  S.l  and  5.2.  Th»  velocity  of  the  yielding  structure  will  ho  ass  weed 
to  be  uniformly  distributed,  that  is,  the  structure  Is  assumed  to  move  data* 
as  a  rigid  bvdy.  This  behavior  correspond?  to  the  stress  field  solutions 
given  In  Section  S.2  in  which  a  region  of  the  grauuler  medium  remains  elastic 
and  moves  essentially  as  a  rigid  S>d*  .  Solutions  for  the  cate  of  a  non-uniform 
motion  cf  the  structure  are  given  in  Section  5.0. 

The  velocity  equation  of  the  detic  system  of  the  second  kinematic 
relation  was  given  in  Section  3.C  in  polar  coordinates,  for  dy/d«  •  tan(0  -  y) 
active,  this  equation  becomes 

sin  ?  ^  -  iSLi  +  -  tin  9*0  (i.4.1) 

or  r  rtf  f 

The  differential  aquation  of  the  character i si  let  which  correspond  to  the  active 
slip  lines  is 


jS  ,  .-.jhll£  (».4>t 

dr  r 

and  the  differential  teUtlon  which  must  be  satisfied  along  the  character  1st Ic* 
is 


dv  «  v  ten  *  (» 


equation  yields 


V 


Q  ten  f 


(4.4.J) 
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where  v^  is  an  arbitrary  constant. 

In  region  0C0  ot  Fig.  5.il,  the  velocity  vector  at  any  point  is 
tangent  to  the  slip  line  given  by  Eq.  (5.4.2);  however,  in  region  000*  the 
velocity  vector  is  vertical.  This  requires  that  a  velocity  discontinuity 
forms  on  the  cha.  acler  i  st  ics  0D  .  Similarly,  the  velocity  is  discontinuous 
across  the  characteristic  0C.  Across  a  discontinuity  of  this  type  only  the 

tangential  components  of  velocity  ere  discontinuous  while  the  normal  compo- 
16 

nerts  are  continuous 

The  velocity  in  region  OCD  just  across  the  discontinuity  00  may  ba 
obtained  in  terms  of  the  velocity  of  the  structure  vq,  by  equating  the  normal 
components  of  the  velocity  vectors  on  each  side  of  00  as  is  shoom  In  Fig. 

5  12 (a)  .  From  Fig.  5.12(a) 


v  ■  v  s  i n (t  ♦  6nn) 
qn  q  00 


Since  the  nv.-aml  components  of  velocity  are  equal, 


or 


V  '  ’»  “’*!  '  'T> 


,  *  V 

¥00 


cot(^  -  *r) 


(5.4.5) 


(5.4.») 


The  arbitrary  constant  in  Eq.  (5.4.'.)  r»y  he  evaluated  at  follow*: 

(t  ♦  •  I)  *•"  9 

vm  •  v.a'’  00  2' 

00  I 


or 


-  y«.(’  r  '  9oo  *  !>  * 


rl  T00 
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Equation  (5.4.4)  now  becomes 


i 


v 


V  e(0+f--r 

W  1 


-  <J0D)  t,n  * 


(5.4.7) 


The  velocity  on  OC  Is  obtained  by  substituting  the  value  of  6  on  OC  Into 
Eq.  (5.4.7). 

•oc  ■  •  s»’  ""  ’  <»•*••> 

The  velocity  in  region  OCC' ,  v  ,  is  obtained  In  teres  of  v_.  by 

-  p  w. 

equating  the  normal  components  of  velocity  across  OC  as  Is  shown  In  Fig. 
5.12(b). 


vp  cos(r  ♦  &oc  ■  ■  y0C  c0*^f  *  (5*4. t) 

Substituting  Cqs.  (5.4.6)  and  (5.4.6)  Into  Cq.  (5.4.0),  one  obtains  the 

velocity  in  region  OCC',  v  as  a  function  of  v  ,  the  velocity  In  region 

P  <l 

000'  . 


v 

P 


sln(r  ♦  0*^) 

vq  .TnTr  4 


e  OC 


®00>  *•"  * 


(5.4.10) 


'•'or  the  stress  field  given  In  Section  5.1,  6gg  •  «/2  and  0^  ■  0. 
Then  Cq.  15.4. 10)  becomes 


V 

p 


4  tan  • 
v^  ton  y  e* 


(5.4  II) 


It  nay  now  be  shown  that  overall  continuity  of  the  velocity  field 
is  satisfied,  feltlplylng  the  velocity  of  OCC'  by  the  length  CC'  gloom  by 
Cq.  (5.1.14),  on<*  obtains 


-  «■  •  % —  •  «•  •  v 


(5.4.11) 


tan  ♦ 


or 


v  •  CC1  »v  •  L 
P  P 

Thus,  a  velocity  field  has  been  found  which  is  compatible  with  the  corresponding 
stress  field. 


6.  RESULTS  AND  CONCLUSIONS  BASED  OH  THE  QUASI-STATIC  THEORY 


6. 1  Comparison  with  Experiments 

The  results  obtained  in  an  analysis  using  the  Geniev  theory  of 
motion  of  a  granular  media  to  predict  the  intensity  of  air- induced  overpressure 
transmitted  to  an  underground  structure  are  shown  in  Fig.  6.1  along  with  experi¬ 
mental  data  obtained  by  Selig,  McKee  and  Vey.  10,22  The  experimental  work  was 
conducted  with  a  dense,  cohesionless  Ottawa  sand  (p  ■  35°)  contained  in  a 
glass-walled  pressure  device.  The  yielding  structure  was  a  4-inch  square  plate. 
It  may  be  readily  seen  from  Figure  6.1  that  the  experimental  pressures  required 
for  fai lure  are  several  times  greater  than  those  predicted  by  the  Geniev 
theory.  The  lack  of  agreement  between  these  two  curves  is  to  be  expected 
because  of  two  major  factors.  For  one,  the  experiments  were  not  actually  per¬ 
formed  under  conditions  of  two-dimensional  plane-strain  as  ass«ed  In  the 
Geniev  theory,  since  the  glass  walls  Of  the  pressure  box  contributed  same 
frictional  resistance  to  slip.  Also,  and  more  significantly,  the  Geniev 
theory  used  here  stiluses  that  plastic  flow  takes  place  at  constant  volume, 
that  is.  the  flow  is  incompressible.  In  the  case  of  granular  materials,  this 
assumption  is  correct  only  a<ter  initial  slip  hat  occurred  and  the  f lorn*  field 
is  established.  ^  For  a  dense  granular  material.  Initial  slip  oust  be 
accompanied  by  an  increase  in  volume  at  the  Individual  particles  "unlock" 
from  their  positions  In  the  dense  state  and  r.da  vrer  *dj  »c*nt  grains  on  the 
slip-plane.  This  intuitively  apparent  behavior  hat  been  proved  conclusively 
bv  experiment  .  As  slip  first  occurs,  there  Is  a  change  In  volume  of 

the  granule-  mats,  and  after  a  certain  strain  the  deformation  settles  down 
to  one  occurring  at  constant  volume.  It  is  In  this  constant  volume  deformation 
that  th.  Geniev  theory  Is  applicable- 
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The  volume  increase  required  for  slip  in  dense  granular  soils  and 
the  lateral  constraint  provided  by  the  soil  mass  itself  hove  a  significant 
effect  on  the  pressure  required  to  produce  first  slip.  Terzaghi  has  Indicated 
that  the  so-called  internal  friction  of  the  soil  is  a  function  of  the  amount 
cf  lateral  expansion  possible  in  the  soil  and  that  the  angle  9  may  vary  between 

25 

wide  limits. 

In  another  paper,  Terzaghi  states  that  for  granular  materials  the 

value  of  9,  the  angle  of  internal  friction,  may  vary  between  the  limits  of 
26 

30  and  54  degrees-  If  in  the  Geniev  analysis  the  value  of  9  is  allowed  to 

vary  in  accordance  with  the  degree  of  lateral  constraint  provided  by  the  soil 

mass  at  various  deptl.s,  a  failure  pressure  vs.  depth  of  burial  curve  Is  obtained 

which  gives  much  improved  agreement  with  the  experimental  results  as  shown  in 

Fig.  6.1.  The  assumed  variation  of  9  with  depth  is  presented  in  Fig.  6.2  as 

the  ratio  of  the  value  of  9  used  in  determining  the  failure  pressure  at  a 

particular  death  to  9  ,  the  value  of  Internal  friction  determined  by  the  usual 
o 

laboratory  tests.  The  value  of  */90  I*  aisumed  to  vary  from  1.0  for  structures 

placed  near  the  grout.*  surface  where  there  >r  a  small  lateral  constraint  to 

1.5  for  a  depth  of  euriel  eque!  to  the  span  length.  At  greatar  depths  the 

value  of  */9  is  atstmwd  to  remain  constant.  Although  the  angle  of  internal 

o 

friction  $  it  tKown  to  v#ry  with  d*pth  in  FI).  1.2,  it  it  not  truttf  at  a 
variable  at  any  on#  particular  depth.  Solutions  Including  the  effect  of 
lateral  constraint  may  be  obtained  by  using  iN.  *j<ee  of  9Q  to  determine  the 
shape  of  the  slip  field  end  then  Incorporating  ths  effect  of  lateral  constraint 
by  determining  the  failure  pressure  based  on  the  value  of  9  given  In  Fig.  6-2 
for  the  particular  depth  associated  with  the  slip  field.  This  approach  it 
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undoubtedly  only  approximate  but  it  does  provide  an  indication  of  the  pressure 
required  for  first  slip  of  the  soil. 

Since  the  pressure  required  to  maintain  slip  in  a  dense  granular 
material  is  less  than  that  required  to  produce  slip,  it  is  clear  that  a 
question  of  stability  is  involved  in  regard  to  the  safety  Of  the  structure 
after  first  slip  h«s  occurred.  Whether  or  not  complete  collapse  of  the 
structure  will  take  place  depends  on  the  intensity  and  timewise  variation  of 
the  loading  pulse  and  the  ductility  of  the  structure.  However,  the  question 
of  stability  does  not  apply  if  the  failure  pressures  are  determined  by  the 
Geniev  analysis  (without  the  lateral  constraint  correction),  since  the  failure 
pressure  is  based  on  that  required  to  maintain  flow  not  s.„,r  it,  and  thereby 
neglects  wh«t  may  be  an  appreciable,  although  highly  undependable,  resistance 
of  the  soi I . 

U  Comparison  with  Previous  Theoretic#!  Studies 

A  theory  of  plastic  flow  of  granular  madia  developed  by  Orucker 
and  Areger  based  on  the  concepts  of  the  plettlc  potential  ties  been 

applied  by  Shield  30  to  problems  of  plestlc  flow  in  granular  soils.  In  this 
theory  plastic  flow  of  a  granular  media  is  airniyi  accompanied  hy  e  volume 
ncreese  (a  property  Oucker  and  Areger  refer  to  es  di'etancy).  However 
in  this  theory  the  voltaae  expansion  continues  et  the  seme  rete  for  all  values 
of  Strain;  a  very  unlikely  event  for  continued  clastic  flow  and  a  behavior 

•  j 

which  experimental  evidence  has  shown  does  not  occur.  Geniev  cites 
experimental  work  on  a  particular  problem  solved  enelyticelly  by  Shield  using 
the  Orucker-Areger  Theory.  The  test  results  indicate  that  the  flow  tra¬ 
jectories  predicted  on  the  bests  of  e  constant  rase  c<  voiiaaetnc  strain  thee  y 
are  not  realized.  Aether,  the  experimentally  determined  flow  field  could  be 
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better  approximated  by  a  constant  volume  (zero  rate  of  volumetric  strain) 
theory  such  as  Geniev's. 

Another  solution  of  the  underground  structure  problem  is  included 
in  a  theory  presented  by  Terzaghi  for  the  design  of  tunnels  at  great  depths.1 
This  theory,  modified  slightly  to  correspond  to  the  specific  problem  considered 
here,  is  based  on  the  assumption  that  the  vertical  stresses  are  uniformly 
distributed  on  horizontal  sections  and  that  the  slip-planes  form  vertically. 
According  to  tills  theory  the  pressure  applied  at  the  ground  surface  to  cause 
fa i  lure  i s  given  by 

p^  •  q  exp  (2  j-  K  tan  qi)  (6.2. 1 ) 

where  K  is  an  empirical  coefficient  which  represents  the  ratio  on  the 

surface  of  sliding.  Terzaghi  recommended  that  the  value  of  K  le  at  least  equal 
to  unity-  Failure  pressure  vs.  depth  curvet  for  Terzaghi's  theory  with 
K  »  I  and  K  •  I.S  are  shown  in  Fig-  6.3  along  with  a  curve  for  the  Geniev 
theory  for  «  «  3S°  and  k  ■  0-  The  experimental  curve  given  in  Fig.  6.1  Is 
also  shown  for  comparison. 

It  may  be  seen  from  Fig.  6.3  that  the  general  trend  of  the  Terzaghi 
curves  does  not  follow  either  the  Geniev  or  the  experimental  curve  except 
at  depth-span  ratios  greater  than  about  I.S,  where  the  Terzaghi  t-nd  Geniev 
curves  predict  very  nearly  the  same  failure  pressure. 

Also  show*  in  Fig.  6.3  is  »  curve  for  a  theory  present*  J  by  ietlg, 
McKee  and  Vey.  the  same  group  who  presented  :’-.r  experimental  results  shot*  In 
Figs.  6.1  and  6.3.  Their  theory  it  based  on  the  formation  of  vor.lcel  slip 
planes  and  a  uniform  distribution  of  vertical  strata  on  horizontal  sections. 

For  the  quasi  - stet  Ic  case  of  e  uniform,  slowly  varying  ove.  .restore,  this 
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theory  is  essentially  the  same  as  Teriaghi's  with  a  linear  variation  of 
failure  pressure  vs.  depth  on  semi- logari thmi c  coordinates  as  shown  in  Fig 
Takagi  has  presented  a  theory  of  granular  soils  deformation  in 
which  the  volume  change  resulting  from  the  plastic  deformation  may  be  taker 
into  account-  Unfortunately,  the  rete  of  volume  change  is  considered  to  bt 
constant  depending  on  the  nature  of  the  deofrmation.  As  a  result,  Takagi  'i 
theory  does  not  account  for  the  real  behavior  of  expansion  or  contraction 
the  start  of  motion  end  subseguent  slip  at  constant  volume. 
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6  3  Li  mi tg t ions  of  the  Theo- 


The  rasu’ts  of  chi s  study  indicate  that  the  Geniev  theory  he* 
several  s . go i f i can*  limitations  when  applied  to  the  problem  of  underground 
structures  in  a  granular  soil  subjected  to  air  overpressure*.  The  fact  that 
the  deformation  is  assumed  to  occur  at  constant  voltme  does  not  allow  the 
volume  change  occurring  at  initial  slip  of  a  dense  granular  material  to  be 
taken  into  account-  This  neglect  of  the  volume  increase  et  initial  slip  is 
particularly  significant  in  the  type  of  problem  considered  here,  as  the 
tendency  fer  a  volume  increase  br  ngs  into  pl«v  *  considerable  '-ateral  con¬ 
strain. nq  effect  ir.  the  soil  which  serves  to  g.  early  increase  the  resistance 
of  the  soil  structure  system  against  slip.  In  contrast  to  this,  a  constant 
voltmw  theory  Should  introduce  little  error  for  problems  in  sdtich  the  soil  i* 
permitted  to  eepaou  at  a  free  surface  such  as  In  the  determination  of  the 
pessive  earth  pressure  of  retaining  -ells  end  the  bearing  capacity  of  footing*. 

Another  weeklies*  of  the  Geniev  theory  common  to  *11  theories  of 
rigid  pii.tir  behavior  is  the  feet  that  no  information  is  evellebte  regarding 
the  formation  of  fe.'u-e  acmes  from  the  pr«vio-«  .lestic  equilibrium  state 
1-  *th*c  word*,  with  rigid-plastic  theories  it  is  not  possible  to  predict  the 
origin  and  trace  the  growth  of  distressed  regions  In  the  medl.»  since  the 
elastic  stress  distribution  prior  to  flow  is  not  available  0"  the  other 
Hend.  with  elesto-plestic  theories  the  formation  of  fa-  Ure  aones  may  be 
premcied  froa  the  elastic  solution  by  incre'-  •»*  txe  loeds  until  the  flow 
criterion  ■*  satisfied. 

An  important  complication  in  theories  of  groovier  media  arise* 
free-  the  fact  that  the  failure  criterion  is  e  function  of  the  m»«n  r-omel 
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stress  a'  a  pcir.t  and  not  siwp'y  »  function  of  the  deviatoric  conponen'  of 
stress  as  is  the  case  with  ideal  plasticity.  For  plasticity  probleos  it  has 
been  shown  that  upper  and  lower  bounds  for  the  fal  lure  load  My  be  obtained 
by  application  of  limit  theorems  32,33,3*'35.  In  addition,  these  theorems 
guarantee  that  any  solution  which  is  both  statically  and  geoeetrlcally 
adei stable  oust  give  the  correct  failure  load.  When  the  yield  criterion 
depends  on  the  Man  noneal  stress,  the  bound  theorems  of  plasticity  do  not 
apply  and  the  failure  load  cannot  easily  be  bounded.  However,  a  theory  such 
as  ten  lev's  which  neglects  a  portion  of  the  soil's  resistance  will  give 
conservative  failure  pressures. 

failure  ei natures  obtained  with  the  teniev  theory  nay  be  Improved 
to  account  appro* ;wat* ly  for  the  effect  of  leterel  constraint  of  the  soil  at 
initial  slip  by  attuning  the  angle  of  Internal  friction  9  to  be  a  functlan  of 
the  .mount  of  sol  I  cover  over  the  roof  of  tho  structure. 

Since  lerge  deformations  of  e  granular  wed  I  a  occur  essentially  at 
constant  voluac,  the  teniev  solutions  should  provide  e  good  est least e  of  the 
resistance  of  the  sol  I -structure  systen  to  sustained  loads  If  the  structure 
has  sufficient  ductility  to  withstand  the  Initial  dlsptecooent. 

The  results  of  this  study  Indicate  that  e  theory  nore  roollotlc 
then  teniev 1 »  Is  needs e  for  tho  analysis  of  the  truly  dynenl «  prob tan  of 
sol  I -structure  Interaction.  This  nore  # labor,  u  -heory  oust  nooo-oorlly 
include  the  effects  of  pan si  on  or  contraction  of  tho  granular  sol  I  as 

notion  takas  p luce  and  the  relation  between  this  velime  change  and  tho 
failure  lew.  in  addition,  the  oode  of  developeent  of  failure  senes  regal  res 
clarification  in  future  work. 
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PIC  2.4  MOHR'S  CIRCLE  OP  9BUUR  RA3S8 


pro.  2.$  Mon's  circle  cp  amssss 
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APPENDIX:  APPROXIMATE  DYNAMIC  THEORY  FOR  SOIL- STRUCTURE  INTERACTION 


A . 1  General 

The  purpose  of  this  Aopendix  is  to  present  a  simple  approximate 
method  of  introducing  time  into  the  Muas i -Stat i c  solutions  obtained  in 
Chapter  S.  in  developing  the  approximate  method  it  is  assumed  that  the  slip 
lines  found  for  the  Muasi-static  solution  do  not  change  appreciably  with 
time.  In  this  way,  a  one- dimensional  treatment  is  possible  which,  for  long 
times,  approaches  the  guasi-static  case  previously  found.  For  simplicity  in 
this  treatment  the  resistance  of  the  structure  was  taken  as  constant. 

For  extremely  snort  times,  of  the  order  of  the  time  for  the  first 
effects  to  reach  the  structure  from  the  surface  above  it,  it  is  n„t  easy  to 
give  an  a  priori  justification  of  the  present  approximate  theory.  The  final 
evaluation  of  this  approximat ior.  will  have  to  await  comparison  with  more 
exact  analytical  solutions  or  with  -xperiments. 

We  consider  here  the  static  problem,  treated  in  Section  S.l,  of  a 
long  plate  below  a  finite  depth  of  soil  layer.  The  so>l  layer  is  supported 
by  a  'igid  plane  except  for  the  pert  carried  by  the  plate.  The  static  and 
the  quasi-static  solutions  of  this  problem  with  a  uniform  static  load  on  the 
top  of  the  soil  layer  is  treated  in  great  detail  in  the  main  body  of  this 
report.  Referring  to  rig.  A. I,  the  velues  of  0  in  the  regions  0A1  and  0C0 
are  0  and  t  ■j,  respectively.  In  the  region  OBD,  B  is  given  by 


B  -  f  -  r; 


<  0  <  (r  ♦  t) 


Equation  of  slip  line  10  is 


1  ,  -  (P  -  r)  tan  9 

r  -y  l  cosec  re 


97 


The  depth  of  the  soil  layer  as  a  function  of  the  span  of  the  plate 


and  the  angle  of  internal  friction  of  the  soil  is  given  by 


h  “  l  (cot  y  +  e"2  tan  ;  y  ■  ^ 


h  =  3  L  tan  «  ;  u-|+| 

o  i  4  2 


A.  2  Basic  Equations  of  Approximate  Qne-0 imensi ona 1  Problem 

The  stress  components  are  expressed  in  the  fol lowing  way: 


"  I  +  sing  cos  2P~]  H(x,t) 


L 


i  -  sin  rp  cos  20 


(A. 2.  la) 


0  »  H(x,t) 

y 


(A. 2.  lb) 


„  «  sin  ;p  sm  20  1  H(x.t)  . 

°xy  L  1  -  sin  9  cos  20J  (A  ‘ 2  ' 

Here  0  ■  0(x , y)  is  known  througnout  the  region  and  H(x,t)  is 

unknown.  We  can  isolate  «  small  layer  of  the  soil  such  as  A^  A^  B(  B^  at 

any  height  x,  and  consider  the  equilibrium  of  that  ylemental  layer.  From 

fig.  A. 2,  it  can  be  seen  that  the  layer  is  acted  on  by  the  vertical  average 

stress  *(x)  M(x , *)  along  A.  6„  and  the  normal  and  tangential  stresses  o  and 

u  u  n 

t  on  the  inclined  faces  A.  A.  and  B.  B_-  The  stresses  on  the  inclined 
n  i  U  I  u 

faces  have  to  satisfy  the  boundary  condition 


t  •  o  tan  » 
n  n 

To  find  the  normal  component  o  ,  we  refer  to  Fig.  A. 3  which  ihowi 

n 

an  elemental  soil  mass  at  the  .  rclined  boundary  with  all  the  stresses  acting 
on  itr  three  faces.  From  F19.  A. 3  we  have 
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X  »  o  4*  +  a  cos  (B  +  y  +  <p)  +o  sin  (B  +  y  +  <p) 

x  ds  xy  ds  x  xy  '  T 


Y  ■  a  4^  +  o  4^*o  sin  (B  +  y  +  <p)  4  a  cos  (B  +  y  4  cp) 
y  ds  xy  ds  y  xy 


a  *  X  cos  (P  4  y  4  <p)  4  Y  Sin  (P  4  y  4  <p) 
o 

2  2 

o„  cos  (5  4  y  +  ?)  4o  sin  (B  +  y  +  q>)  +  a  sin  2(p  4  y  4  ip) 
y  xy 


Using  Eqs-  (A. 2.1)  and  simplifying,  wc  obtain  the  normal  and 
tangential  stresses  to  be 

M(x,t) 


( A.  2 . 2a) 


( A.  2 . 2b) 


The  total  vertical  component  c  testes  is 


He  re , 


f  +  4  *<*) 

*  (x)  •  *(*)  /  ,  T  O)  dy 

7  i(x) 

whe  re 

r  ff»  .  +  »!»  T  HLi|1 

I  I  -  sm  ®  cos  2pj 

/(*)  •  length  of  tne  element*!  layer  at  a  height  x. 

X  *■  body  force  intensity 

a  *  acce  le  rat  i  or.  of  the  layer, 

* 

*  t-  ♦  v  r~  .  since  x  is  an  tulerian  coordinate, 

•I  lx 

v(*,l)  •  average  velocity  of  the  layer 


Simplifying  the  above  expression,  we  obtain  the  equation  of  motion 
of  the  layer  >n  the  form 


X 


(x) 


1 

J 


>v 


(A.  2.3! 


where 


the  velocity  v(x,t)  has  to  satisfy  the  continuity  equation  which  is  0>  the 
for* 


<*■*■*> 

Here  we  have  two  u.vt>i*®ens,  K(x,t)  and  v(x,t),  and  two  equations,  (A.2.J)  and 
(A. !  41  to  obtain  the*- 


too 


To  obtain  the  function  \|r(x)  ,  we  have  to  evaluate  the  integral 


♦  (x) 


"777 


7“/_  I 


+  ^  i(x) 
i(x) 


r  (0)  dy 


This  is  difficult  to  do  exactly  since  0(x,y)  is  a  complicated  function  of 
x  and  y  and,  also,  the  limit  y  «  t  i(x)  is  actually  a  transcendental  function. 
Instead,  we  can  find  an  approximate  expression  for  the  average  value  of 
0(x,y)  over  the  range  *h  <  x  <  0  as  a  function  of  x  only.  Then, 

0(x,y)  ■  0(x)  and  r(p)  -  r(P) 

We  can  obtain  the  function  t(x)  easily  as  follows: 


*(x) 


+  j  iU) 


1W  dv  * r(5) 


Now,  P(x)  has  to  satisfy  the  following  conditions: 


P  -  0. 

f  -  0  x  -  0 

(A.S.Ia) 

p.f 

X 

1 

• 

* 

o 

■ 

(A. 3  lb) 

P  .  V. 

x  ■  -h 

• 

(A.  3.1c) 

We  can  take  0(x)  as  a 

obtain  the  constants  to  satisfy 

polynomial  in  x  with  fi 

tht  conditions  (A»3«l). 

ve  constants  and 

Then  we  oet  the 

following  expression  for  P(x) : 


where 


Now, 


(A. 3. 2) 


foU)  dr(&) 
dx  IT  • 


m  .  f  *  >iB  9  >ln  ^  1 

1(1  -  tin  v  tin  2b)  J 

We  need  thit  derivative  of  f(x)  in  evaluating  the  function  ((x) . 

The  function  j(x)  which  givct  the  width  of  the  elewental  layer  at 
any  height  x  it  required  to  evaluate  the  function  ((x)  which  It  a  variable 
coefficient  of  Eo-  (A. 2. 3).  In  cartetian  coordinate!,  If  we  eatign  any  value 
for  x,  we  have  to  coafUe  i(x)  at  a  tolution  of  a  trantcendentel  equation,  at 
the  equation  of  the  curve  10  in  Fig.  A. I  it  known  In  polar  coordinate!. 
Inttead,  we  can  ettignveluet  of  angle  0  In  the  range  y  and  (y  +  -j)  and  then 
compute  the  valuet  of  J(x)  in  a  very  tlnple  nh'nnvr  from  the  following 
equation!; 

1(10  ■  Tir;  [""  *  ’ r> “"  *  ] 

•  -  t  +  r  (A.3  3) 

r  <  P  <  y  •>  -8 
-  -  t 
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Tnese  equations  are  useful  to  compute  the  functions  y(x)  ana  £vx)  as  functions 
of  9. 

A.  4  Boundary  and  Initial  Conditions 

On  the  top  of  the  soil  layer  a  uniform  load  which  is  a  function  of 
time  acts  vertically  down.  At  the  bottom  of  the  soil  layer,  the  component 
of  the  stress  in  the  soil  normal  tc  the  plate  at  the  interface  will  be  the 
load  on  the  plate  so  that  this  is  resisted  by  the  inertia  of  the  plate  as 
well  as  its  yield  resistance.  It  is  assumed  that  there  is  no  separation 
between  soil  and  the  plate  along  the  interface  so  that  the  velocity  of  a 
point  of  the  plate  is  the  same  as  the  adjacent  velocity  of  the  soil.  Hence, 
the  boundary  conditions  may  be  defined  as  follows: 
on  the  top  where  x  »  -h 

*  P(t) 

at  the  bottom  where  x  »  0 

°x  *  °p  ~V'^5t~^  + 

where 

Op  ■  density  of  the  material  of  the  plate, 

h_  *  thickness  of  the  plate 
P 

K  ■  yield  resistance  of  the  pi*'.-. 

Using  Eqs.  (A. 2.1)  the  boundary  conditions  may  be  written  at 

H  (-h.t)  -  (}  !  ffi"*)  P(0  *  -  -h  (A. 4.  la) 
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•<*  * 


x  “  0 


(A-4.lb) 


~  r 

s' 


Tne  initial  conditions  may  be 


H  (x,o)  ■  H  (x)  and  v  (x,0)  •  v  (x)  at  t  •  0 


(A. 4. 2) 


A- 5  Method  of  Solution 


Equations  (A. 2. 3)  and  (A. 2. 4)  are  quasi  linear  with  two  unknowns 
H(x,t)  and  v(x,t)  and  they  can  be  solved  by  the  method  of  characteristics. 
The  density  of  the  soil  p  is  taken  to  be  a  function  of  pressure,  so  that 


»'  »  -g 


Substituting  these  expressions  into  Eq-  (A. 7.4)  we  obtain 


(A.S.l) 


ciH  dH  dH  d< 

St  "St  ’&7t 

gy  ,  dv  dx 
dt  "  dt  *  dt 

Substituting  these  aviations  into  (A. 2. 3)  and  (A.S.l)  and  simplifying,  we 
obtain 
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(A. 5. 2a) 


.  dH  .  f  dx  1  dv 

*  ^  +  p  i v  *  dt  J  S 


dv  1 

ox  -  CH  -  p  df  j 


r  dxl  dH  .  o(H)  dv  _ 
dt  §x  p'  (H)  S< 


djfl 

dt 


(A- 5. 2b) 


From  Eqs.  (A. 5. 2)  we  can  solve  for  the  partial  derivatives  of  H(x,t)  and 
v(x.t)  to  get  the  characteristics  and  the  equations  to  be  integrated  along 
them  as  fo!  lews: 


IS" 


-  tH  *  P  £) 

»<«  - 

(.  HO) 

v  dt' 

o(H) 

P(v  -  g) 

*(x) 

o(H) 

pW 

(A. 5. 3) 


By  setting  the  denominator  equal  to  zero  and  slwplfying,  we  find 
the  characteristic  equations 


dx 

dt 


v  t  u(H) 


(A.S.4) 


where 


u(h)  -  \JtU-;  ^ 

By  setting  the  numerator  equal  to  zero,  we  find  the  equations  to 
be  integrated  along  the  characteristics: 

+  ;(*)  H  -  p  xj  -  0  (A-5-5) 


♦  (*)  $  * 
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Hence,  we  have  the  following  set  of  ordinary  differential  equations: 


dx  - 


u(h; 


v  +  u(H)  j  dt  -  0 

dH  +  o(H)  dv  +  ;(x)  H  -  pxj  dt  -  0 


>  (A. 5. 6a) 


dx  -  £v  -  n(H)j  dt  ■  0 


j-j-  dH  -  p(H)  dv  -  {(x)  H  -  pxj  dt  -  0 


>  (A. 5.66) 


We  have  thus  reduced  the  quasi  I  inear  partial  differential  equations,  (A-2-3) 
and  (A. 2. 4)  to  a  set  of  ordinary  differential  equations  which  can  6«  solved 
numerical ly  by  the  method  of  finite  differences- 


A. 6  Method  of  Finite  Differences 

Equations  (A. 5. 6)  may  be  written  in  f Inite-di fference  form  for 
the  purpose  of  numerical  computation.  The  first  step  Is  to  transform  iqs. 
(A. 5. 6)  to  curvilinear  coordinates  of  the  characteristics.  Let  the  first 
and  the  second  character! st 1 cs,  respectively,  be  given  by 

U  »  U(*,t)  ■  constant 
Z  •  Z{x,t)  »  constant 

We  can  take  these  as  curvilinear  coordinates  and  transform  fqs.  (A. 1.6) 
from  xt-plan*  to  U2-plane  by  treating  a,  t,  H  and  v  as  functions  of  0  and  I. 
Along  the  first  system  of  characteristics  we  have 

dx  -  ^  dZ  ;  dt  - 

«»-$“!  «»  -  |?  « 
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Similar  relations  are  written  for  the  U-characteristics. 

With  the  help  of  these  relations,  Eqs-  (a-5.6)  may  be  written  in  the  UZ-plane 
in  the  form: 


Sx 

dZ 


> 


C  H  -  DX 


0 


(A.6.U) 


to  *  [v  -  to  ■  0 


^hT  $>  ■  e  to  “  [t  H  '  px]  $5  "  0 


(A- 6.  lb) 


Now,  the  derivatives  can  be  written  at  differences  over  a  finite  length  as 
f ol lows: 


£ 


V<  ‘.V»i 

£*J 


f 


Ji xL 


as  shown 


following  this  rule,  Eqs  (A.S.6)  are  written  at  difference  equations 
below: 


(\.i  *  Vn1 

*  [\,t- 1  +w{Mk,d- 

r»]  (tk.i  *  ik,j-»)  * 0 

(A, 6. 2a) 

<YI  * 

*  [Vi.!  ’  u(Hh-i. 

l>]  {tk,4*  «h 

(A.«2b) 

Yl-I 

“(Hh,J-l)  k,< 

‘  Vi- I*  *  °k,i-l 

(vk,l*  Vk,4.l)  4 

(A.».2c) 

k ,  i*  I 


°k,|.|  *k,J.| 
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(A. 6. 2d) 


JTThJ^  i)  (Kk,i  '  Hk- 1,1^  *  pk-  I , i  (vk,i  ‘  Vl,J) 

?k-l,l  *k-l,4  ‘  pk-  I ,  i  *k-l,i]  At2  *  0 

Equations  (A. 6. 2a)  and  (A. 6. 2b)  give  the  values  of  x  and  t  at  any 
point  (k,i)  in  the  l>Z-plane  with  the  known  values  of  x,  H  and  v  at  the 

previous  points  (k-l,i)  and  (k,4-l).  Similarly,  F.qs.  (A. 6. 2c)  and  (A.6.2d) 
give  the  values  of  H  and  v  at  the  point  (k,i)  with  the  known  values  at  the 
previous  points  (k-l,i)  and  (k,J-l).  Using  these  relations  a  computational 
scheme  may  be  easily  arrived  at. 

The  network  on  the  U2-p!ane  is  shown  In  Fig.  A. a.  pg  represents  the 
initial  conditions-  In  the  region  A8C,  along  the  line  Al,  we  know  the  values 
of  x,  t,  H  and  v.  the  given  initial  conditions.  At  points  along  t  ■  0,  we 
have  the  following  data: 


Aoints 

(6.1) 

(S  Z) 

(♦.3) 

(3.4) 

(2.S) 

0.6) 

t 

0 

0 

0 

0 

0 

0 

«/h 

0 

-0.2 

-0.4 

-06 

-0.8 

-1.0 

H 

*61 

HS2 

*43 

*34 

*25 

*.* 

V 

W6I 

¥S2 

*43 

v34 

V2S 

’18 

We  may 

compute  x,  t. 

M  and  v  i 

■t  various  points 

in  the  region  AtC 

in  the  following  manner; 
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(6.2)  (5,3)  (4,4)  (3,5)  (2,6) 

(6,3)  (5,4)  (4,5)  (3,6) 


where 


(6,4)  (5,5)  (4,6) 

(6,5)  (5.6) 

(6,6) 


At  each  of  these  points  we  have  the  following  four  equations: 

(*66  -  X6S}  *  [V6S  +  *(IW]  (t66  ’  W  "  0  fA63al 

*x66  ‘  *56^  ‘  [v56  *  ^H56*]  **66  *  *56*  “  0  (A. 6. 3b) 

JTO  (M66  *  M65}  +  °65  (w66  *  v65}  + 

“  (A.  6-3c) 

*^65  “65  '  p65  W  •  A*1  *  0 


♦56 

uTh^T 


(H, 


66 


V 


°56  (v66  ‘  y56*  * 

(A3. 3d) 

^56  *56  *  *S6^  '  •  ® 


61 1  “  (t66  *  *65^ 


At2  "  (,66  ' 


ir.  II.*  region  ACO,  we  My  proceed  along  the  points  shown  bel»: 
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(2,7) 


(3.7) 


(4.7) 


(5.7) 


(5.7) 


(3.8)  (4,8)  (5.8)  (6,8) 

(4,9)  (5,9)  (6,9) 

(5,  10)  (6,10) 

(5.11) 


To  compute  the  '-alues  a t  the  point  (2,7) : 

(x27  *  X26)  ‘  [ v26  4  (t27  '  t26>  "  0  (A#>4*) 

*j7  ■  -h  (A. 6. 4b) 

{H27  •  M2«)  +  °26  {v27  *  V2*}  4  (A*4c) 

4  «2*  M28  *  p2*  V  (t27  *  *2.)  *  0 
H27  «  ‘><t27>  *•••*> 


Thete  equation*  ere  to  be  uted  at  points  (2,7),  (3.8),  (4,8), 

(5,10)  and  (6,11),  which  are  the  boundary  point*.  At  the  other  point*,  the 
vaM  equation*  a*  In  the  region  AIC  are  to  be  used. 

In  the  region  ICE,  the  value*  at  point*  are  computed  In  the  foilwr¬ 


ing  manner: 


(7.3) 

(7.4) 

(7.5) 

(*.4) 

(*.5) 

(9,4) 

(9.5) 

(9.*) 

(10,5)  (10.*) 

(»1. *1 


(*.*) 


(7.4) 


no 


1 


To  compute  the  values  at  the  point  (7,2)  we  use  the  following 

equat ions: 


^x72  "  X62)  ‘  (^62  *  m^H62* 


x  *0 

*n  u 


(t72  *  t62)  *  0 


(A. 6. 4a) 
(A. 6. 4b) 


*f2 


u(H,,)  ^H72  *  n62;  '  M62  v’72  "  ’62' 


(A -6.4c) 


62' 


*  U62  H62  “  p62  X62}  *l72  *  t62)  "  0 

rV,„  v. 


"7!  ■  [»r  fcp  *  ']  (*'6  4',) 


These  equations  should  be  used  at  points  (7,2)  ,  (8,3)  ,  (9,4)  , 

(10,5)  and  (11,6)  as  the  boundary  points.  At  other  points,  the  same  equations 
as  in  region  A8C  are  used. 

In  the  region  OCEF,  the  unknowns  at  various  points  are  caagsuted  as 

fot lows 


(7.7) 

(8.7) 

(9.7) 

(10,7) 

(U.7) 

(?.8) 

(8.8) 

(9.8) 

00.8) 

(H.8) 

(7,9) 

(8.9) 

(9.9) 

(10,9) 

(11.9) 

(7,10) 

(8.10) 

(9.10) 

(10.  .0) 

(n.io) 

(7.11) 

(8.11) 

(9.11) 

(10.11) 

(11.11) 

At  all  these  eoints,  sane  equation*  as  In  the  region  A8C  are  used. 
The  reg>  S"  0*6  ■»  similar  to  the  region  ACC  and  the  region  tfl  is 


similar  to  the  region  IC(  and  to  on. 


I 


*(*)  HU,t) 


mb.  a.* 


id 
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